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Abstract 



A total and differential cross-section calculation for the Electroweak Standard Model re- 
action e'^e" — >• (Z^Z^) — > /i/i/2/2(7) including the effects of the finite Z^ width and, 
for the first time, all QED initial state radiative corrections (ISR) is presented. A semi- 
analytical technique is used: All angular phase space integrations, five for the tree level 
process, seven if ISR is included, are carried out analytically. The remaining phase space 
variables arc the invariant masses of the two decaying Z^ bosons for tree level plus the 
reduced center of mass energy squared in the ISR case. Invariant masses are submitted to 
high precision numerical integration. Semi-analytical and numerical results for e+e" cen- 
ter of mass energies between ^/s — 130 GeV and 1 TeV are reported. It is shown that 
the radiatively corrected cross-section splits into a universal, leading logarithm part with 
the tree level cross-section factorizing on one hand and a non-universal contribution on 
the other hand. A generalization to the case of Standard Model neutral boson pair pro- 
duction, e+e" — > Z°7, 77) — > /i/i/2/2(7), is worked out in detail. The relevance 
of the presented calculation for the phenomenology of e'^e~ collider physics is discussed, 
especially for W and Higgs physics both of which will be important topics in the near 
future. All significant technical details of the calculation are explained. 



Zusammenfassung 

Im Rahmen des Standardmodells elektroschwacher Wechselwirkungen wird eine Berech- 
nung totaler und different ieller Wirkungsquerschnitte fiir die Reaktion e~^e~ (Z^Z^) 
/i/i/2/2(7) vorgestellt. Dabei wcrdcn die Effekte der endlichen Bosonlebensdauer und 
erstmalig die voUstandigen Eingangszustands-QED-Strahlungskorrekturen (ISR) beriick- 
sichtigt. Es wird eine semianalytische Technik angewandt: Samthche Winkelvariablen des 
Phasenraums, fiinf fiir die Born- Approximation, sieben im ISR-Fall, werden analytisch in- 
tegricrt. Die verbleibenden Phasenraumvariablen sind invariante Massen, namlich die in- 
varianten Massen der zcrfallcnden Z^'-Bosonen und im ISR-Fall zusatzlich die reduzierte 
Schwerpunktsenergie. Uber die invarianten Massen wird mit einem sehr prazisen nu- 
merischen Algorithmus integriert. Es wird iiber semianalytische und numerische Ergeb- 
nisse fiir e+e"-Schwerpunktsenergien zwischen y^^lSOGeV und 1 Te V berichtet. Dabei 
wird gezeigt, dafi der Wirkungsquerschnitt mit ISR in einen universellcn Anteil fiihrender 
Logarithmen mit dem Born- Wirkungsquerschnitt als Faktor einerseits und einen nichtuni- 
versellen Beitrag andererseits zerlegt werden kann. Eine Verallgemeinerung der Resultate 
auf die Standardmodell-Paarproduktion neutraler Bosonen, e'^e~ — > Z°7, 77) — > 

/i A 72/2 (7) wird detailliert ausgearbeitet. Die phanomenologische Bedeutung der vorge- 
stellten Rechnungen fiir die e~^e~ Beschleunigerphysik wird diskuticrt, insbesondere die 
Bedeutung fiir die in naher Zukunft wichtigen Themen der W- und der Higgs-Physik. 
AUe wesentlichen technischen Einzelheiten der Rechnung sind erlautert. 
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Chapter 1 
Introduction 



It is the goal of modern physics to achieve reduction of nature's complexity to as few as 
possible fundamental principles. The method is increasing abstraction, the tool mathema- 
tics, the result unification of, up to that date, distinct phenomena. An important success 
pointing in that direction was Maxwell's formulation of the theory of electromagnetism in 
1864 The electromagnetic theory was not only the origin of the theory of relativity f^, 
but it could also be cast into a Lorentz-covariant form. In the late twenties, originating 
from Dirac's work [0], electromagnetism continued its career as a quantum field theory 0] 
which was shown to be U(l) gauge invariant by Weyl 0. Quantum Electrodynamics 
(QED) was born, when Schwinger reformulated the electromagnetic quantum field theory 
in a Lorentz-covariant way . Shortly thereafter the well-known Feynman diagrams were 
introduced 0, and it only took another year to establish the mathematical foundations 
of Feynman's theory |Q. Today, QED is one of the best-tested theories we know. A 
manifestation of this success is the prediction of the anomalous magnetic moment a^. of 
the electron up to 5a(./ae = 2.4 x 10^® from perturbative QED 0] in agreement with 



experiment |T^. Other QED predictions matching high precision experiments are the 
anomalous magnetic moment of the muon [11| and the Lamb shift of the hydrogen 



atom [y, |T2[. This list could be continued. The successes of QED are clear evidence for 
the quantum corrections as obtained from relativistic quantum field theory. Perturbative 
QED is a very important ingredient of the Standard Model of Electroweak Interactions 
(SM). 

The second main ingredient of the Standard Model are the weak interactions which en- 
tered history in 1896 when Becquerel discovered radioactivity in the form of nuclear 
/9-decay |]l3l- In 1934 Fermi developed a "theory of weak interactions" for /5-decay which 
was constructed along the lines of QED WM. The structure of weak interactions was fur- 



ther clarified with the theoretical consideration of parity violation by Lee and Yang |15 



Soon thereafter, parity violation was experimentally confirmed in the famous Wu ex- 



periment [TO. This was when Fermi's original pure vector currents were replaced by a 



(V-A) structure which is maintained until today Jl^. Although very successful, Fermi's 
theory of weak interactions is not renormalizable and runs into a unitarity problem at 
center of mass energies around 600 GeV. This problem was cured by the introduction of 
gauge bosons that mediate the weak interaction. The short range of weak interactions 
requires massive gauge bosons via the Yukawa theory [0 . However, for the sake of gauge 
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symmetry, bosons are not allowed to have explicit mass terms within the framework of 
a Yang-Mills field theory [0. A remedy to this mass prohibition was found with the 
Higgs(-Kibble) mechanism of spontaneous symmetry breaking which "generates" gauge 
boson masses that are proportional to the non-vanishing vacuum expectation value of a 
newly introduced scalar field [^. In addition, the Higgs mechanism enables the gauge 
invariant introduction of chiral fermion masses into a Yang-Mills theory. In 1971, it could 
be proven that Yang-Mills theories with gauge boson masses generated by spontaneous 
symmetry breaking are renormalizable pT| . 

It has been a milestone of modern physics when, in parallel to the development of the 
Higgs mechanism in the sixties, a renormalizable, non-abelian SU{2) l x U{1)y gauge field 
theory for the unified description of weak interactions and QED could be found |^ which 
is today called the Standard Model of Electroweak Interactions. The groups SU{2)l and 
f/(l)y are associated with the weak isospin and the weak hypercharge symmetries. Via 
the Higgs mechanism three of the four gauge bosons of the model, namely the and the 
are endowed with masses, whereas the photon as the mediator of the electromagnetic 
force stays massless as required by the infinite range of QED forces. In the framework of 
the Standard Model, the Higgs mechanism entails the introduction of one scalar boson, 
the so-called Higgs boson. The photon has a pure vector coupling to electromagnetically 
charged particles, bosons have maximally parity violating (V-A) coupling struc- 
ture, and the Z° bosons comprise both (V-A) and (V-l-A) couplings. In the 1950s and 
1960s many weak decays of hadrons were experimentally studied. With respect to these 
observations it is a crucial feature of the Standard Model that it is able to naturally in- 
corporate the hadronic constituents, the quarks p3|, via the Glashow-Iliopoulos-Maiani 



mechanism ||2^. This incorporation was based on the idea of a unitary transformation 
of the quark flavor mass eigenstates into weak eigenstates which is today known as the 
Cabibbo-Kobayashi-Maskawa matrix 03. The present version of the Standard Model is 



made anomaly- free by the introduction of color p6|, a concept originating from Quantum 
Chromodynamics 



The experimental confirmation of the Standard Model started in 1973 with the discovery 
of the predicted neutral current interactions in the Gargamelle bubble chamber at the 
CERN proton synchrotron, where i^j scattering on electrons and nucleons was investi- 



gated Soon the Standard Model was further supported by the discovery of the weak 
isospin partner of the strange quark, which was called "charm" [^. The breakthrough 
came with the direct observation of the gauge bosons and Z*^ by the UAl and UA2 
experiments at the CERN SppS [0. With the advent of the SLC and LEP e'^e~ colliders 
for resonance production of Z^ bosons in e^e~ annihilation, the Standard Model entered 
its precision test era. From LEP and SLC, the mass Mz and the width Tz of the Z^ as 
well as the Weinberg mixing angle sin^^g^^ are very accurately known [pl[] : 



Mz = (91.1887 ± 0.0022) GeV 
Tz = ( 2.4971 ± 0.0032) GeV 
sm^e[jf = 0.2315 ±0.0004 . 

These impressive numbers, however, are only the peak of an iceberg of electroweak pre- 
cision measurements all of which agree with the Standard Model. The measurements at 
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LEP have reached a level of precision which enables to test radiative corrections not only 
in the QED sector of the theory. Also the truly electroweak radiative corrections and thus 
the quantum corrections of the Standard Model are now established |^. The renormal- 
izability of the Standard Model ensures the sensibility of the perturbation expansion of 
these quantum effects. The latest success in the confirmation of the Standard Model is 
the discovery of the last missing fermion of the Model, the top quark, by the CDF and 
DO collaborations at the Fermilab "Tevatron" pp collider ^^f\. It is very noteworthy that 
the direct measurements of the top quark mass are in good agreement with indirect mass 
determinations obtained from precision experiments at the pole ||3l|, |3^. Precision 



measurements can also severely constrain new physics beyond the Standard Model, and it 
is discussed, whether they can put limits on the Higgs mass, if the precision of top quark 
mass measurements is further improved. 

An upgraded version of LEP, commonly referred to as LEP2, will pass the threshold for 
the pair production of heavy gauge bosons in the very near future and operate at center 



of mass energies between 175 and 206 GeV This opens the opportunity to study the 



properties of the W boson in W pair and single W production processes. Furthermore, 
the non-abelian coupling structure of the electroweak theory will be investigated directly 
for the first time through the observation of trilinear '~fW~^W~ and Z^W^W~ interac- 
tions. Recently, the different non-standard trilinear boson couplings have been subject of 
several studies If LEP2 does not yield the discovery of the Higgs boson, it will at 



least increase the lower Higgs mass limit to approximately the machine energy minus 100 
GeV, i.e. \fs — 100 GeV [^, 0. The above-mentioned physics tasks at LEP2 involve the 



analysis of four-fermion final states. At present, the latter are heavily studied, and several 
experimental analyses were already carried out at LEPl I^U], This thesis is a contri- 
bution to the clarification of the physics of four-fermion final states in e^e~ annihilation. 
A study of the process 



e^e- ^{Z'Z')^ hhhh{^) , h^e^,ue (1.1) 

will be presented, including the first complete treatment of initial state radiation in first 
order perturbation theory [JT|. Process (1.1) is not only observable at LEP2, but also 



interesting at planned higher energy e^e~ colliders. Process (|1 . 1| ) yields a higher order 
test of the Standard Model and can easily be extended to nonstandard neutral current 
physics with additional neutral bosons replacing the Z°'s. Another strong motivation for 
this study is the importance of process ( |1 . 1| ) for W pair physics and for Higgs searches at 



LEP2. The relevance to W pair physics is twofold. Process ( |1.1|) does not only consti- 
tute a double resonance background to some decay channels in W pair production, but it 
was also important to confirm the generality of properties found for W pair production 



in the framework of the "current splitting technique" [^2|, As far as Higgs searches 



at LEP2 are concerned the Z° pair double resonance background to four-fermion final 

'''Since the determination of the number of hght neutrino generations as N^, = 2.987 ±0.016 there 
is reason to beheve that only three generations of quarks and leptons exist. Thus, in the framework of the 
Standard Model, the observation of the top quark together with the indirect evidence for the r neutrino 
justifies the statement that the last missing fermion has been discovered. 

■^It should be mentioned that non-standard trilinear boson couplings are heavily constrained by uni- 
tarity requirements |^. 
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states with a b quark pair must be taken into account. This thesis presents sufficiently 
precise results for studies of the Z° pair double resonance background to Higgs searches. 
Further importance of the presented results on process ( |1.1| ) lies in their straightforward 
generalizability to include new heavy neutral bosons. It should be investigated in how 
far Z'^Z' pair production in e"'"e~ annihilation can put competitive limits on the mass 
and/or couplings of a new boson Z'. On-shell Z^ pair production has been discussed 
long ago [0. Numerical calculations including all 0{a) electroweak corrections^] except 
hard photon bremsstrahlung were reported in for on-shell and in for off-shell 



Z° bosons. In this thesis, initial state radiation (ISR) to Z^ pair production will be 
treated completely, including hard bremsstrahlung corrections. Contrary to the more 
common approaches to process ( |1 . 1|) and other neutral current processes with Monte 
Carlo phase space integration this thesis follows a "semi-analytical" approach 

where all angular phase space variables are integrated analytically and the remaining 
invariant masses are subjected to a high precision numerical integration. Thus the preci- 
sion of the cross-section calculation is physics dominated, meaning that the error on the 
computed cross-section essentially originates from theoretical systematic errors such as 
the disregard of other radiative corrections or the small final state mass approximation. 
Another advantage of semi-analytical calculations is their transparency, i.e. the fact that, 
at tree level, they yield very compact expressions for cross-sections and invariant mass 
distributions. Even with initial state radiation included, semi-analytical cross-section for- 
mulae allow to identify physically interesting properties. Having in mind that, using the 
semi-analytical method, it is difficult to implement experimental cuts on all relevant angu- 
lar variables it is clear that semi-analytical and Monte Carlo calculations of four-fermion 
final states complement each other. However, it should be born in mind that experimen- 
tal cuts on invariant masses and the boson scattering angle can be easily implemented in 
semi-analytical calculations, and it is also possible to obtain semi-analytical distributions 



in many interesting variables |50 



Despite its stunning successes, the Standard Model has conceptual drawbacks, and it is ex- 
pected that there will be new physics at higher energies. Firstly, and more philosophically, 
the Standard Model relies on "too many" arbitrary assumptions and parameters. This 
fostered the hope to find an underlying theory that can e.g. explain the origin of mass, 
of electric charge, or color. Secondly, one would like to find a quantum field theory that 
unifies the electroweak with the other known interactions, namely the strong interaction 
described by Quantum Chromodynamics (QCD) and gravitation. Thirdly, the fact that 
the Standard Model is not asymptotically free suggests that it is a kind of "low energy" 
effective theory of a more fundamental one In addition it is known that, within the 



framework of the Standard Model, the scattering of longitudinally polarized W bosons 
violates the Froissart limit for the proper high energy behavior of cross-sections. Finally, 
there is the well-known fine-tuning or naturalness problem of the Standard Model which 
is related to a quadratic divergence in the radiative corrections to the Higgs mass. This 
problem necessitates new phenomena at a scale of 0(1 TeV) ||5ll|. A candidate to solve 
the fine-tuning problem is Supersymmetry which introduces a fundamental symmetry be- 
tween bosons and fermions and thus, via the relative minus sign in fermion loops compared 
to boson loops and via supersymmetric relations between masses and couplings, cancels 



^That is the electroweak corrections of first order in the fine structure constant a = / (47r) 
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the problematic quadratic divergence. Supersymmetric theories are free of quadratic di- 
vergences |5l|, A further interesting class of theories beyond the Standard Model are 
String Theories. In spite of their mathematical difficulties, String Theories are attractive, 
because they are finite to all orders of perturbation theory and ultraviolet divergence 
free, because they rely on few free parameters, because they may be able to explain mass 
hierarchies, and because they yield coupling constant unification. In addition. String The- 
ories comprise gravitation, possibly solve the problem of the non-vanishing cosmological 
constant [Q, and they may also be on their way to reach phenomenological predicta- 
bility [|^. As a final remark it must, however, be stressed that so far no experimental 
evidence for physics beyond the Standard Model exists. 

This thesis is organized as follows. In Chapter ^ the tree level results for process ( |1 . 1| ) 
will be discussed, followed by the presentation of QED initial state radiation (ISR) results 
in chapter |^. In chapter ^, results from preceding chapters are generalized to include s- 
channel exchange photons. Following this body of the thesis a set of appendices contains 
notations, definitions, and technical details of the presented calculations. Appendix ^ 
gives the metric for this thesis and appendix introduces the Standard Model Lagrangian, 
the Feynman rules, and the renormalization that were used. Appendix |C| contains the 
definitions of polylogarithms and in appendix the 2— and 2^5 particle phase space 
parametrizations needed for the evaluation of cross-sections at tree level and with ISR 
are given. In appendix 0, all important details of the calculation of the tree level and 
ISR cross-sections are collected. Appendix |F| presents techniques for the evaluation of 
loop integrals as encountered in the virtual initial state QED corrections to process ( |1.1|) 
and appendix ^ contains all analytical integrals required for the angular phase space 
integrations to be carried out in the cross-section calculation for process ( |1.1| ). 
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Chapter 2 

The Pair Born C ross-Section 



At Born or tree level and for a center of mass energy ^/i, process ( |1.1| ) is described by the 
two Feynman diagrams depicted in figure |2.1| . The corresponding cross-section is derived 
in appendix |E.1| . It can be written as a simple double convolution formula 

o'^is) = J dsi2Pz{su) J dsuPzisu) x 2 ■ BR{1) ■ BR{2) x 

X (xf (s;Si2,S34) (2.1) 

with branching ratios BR{1) and i?i?(2)|^ of the boson into the fermion-antifermion 
pairs /i/i and 72/2 respectively, with invariant fermion pair masses S12 and S34, and with 
the Z^ boson Breit-Wigner densities 



pz{Si 



Sij Tzisi 



IT Si 



S{s 



Ml) 



(2.2) 



ijj 



In the limit of on-shell Z^ bosons their widths vanish and the Breit-Wigner densities are 
replaced by Dirac S distributions. The Z° width is given by 



(2.3) 



where Mz represents the Z° mass, is the Fermi coupling constant, and gjj^^ and (7^^ are 
the vector and axial vector couplings of a fermion-antifermion pair // to the Z^ boson 
(see appendix [B.2| ). The color factor Nc{f) is unity for leptons and three for quarks. 
Using left- and right-handed couplings as introduced in equation (|E.3|) 



^e+e-Z 



R. 



e+e-^Z 



4 sin6'„, cos6'„ 



4 sin6'^ cos6'^ 



■0 



2 sin^e,, 



■ 2 sin^^,, 




^The factor 2 conventionally introduced with the branching ratios, simplifies comparisons with on-shell 
pair production. This is seen from equation (2.7) and was explained in some detail in reference |41 
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Figure 2.1: The Born level Feynman diagrams for off-shell pair production, the so- 
called 'Conversion' (or 'Crab') diagrams. Left: t-channel. Right: u-channel. 



o"4 (s; Si2, S34) is obtained after fivefold analytical integration over the angular degrees of 



freedom of the 2^4 particle phase space (see appendix p.l|) 



erf (s; S12, S34) 



(Ll + Rf)-G't\s-s^,,s,,) 



(2.4) 



with the sub-index 4 indicating that the underlying matrix element contains four resonant 
Z° propagators. In af , final state fermion masses are neglected. The factor in front of ^4"*"" 
comprises phase space factors and couplings to the initial state current. The kinematical 
function ^4'''"(s; S12, S34) for the summed contributions from the t-channel, the u-channel, 
and the t-u interference is given by 



S] S12, S34J 



[Sl2 + S34) 



C 



5l2 ~ •S34 



B 



(2.5) 



with 



A 



C{s] S12, S34) 

A(s,Si2,S34) 



In — 

s 



S12 - S34 



S12 - S34 



^4 

ence 



A(a, 6, c) = + + - 2ab - 2ac - 2bc . 

was rederived and fully agrees with earlier results M, E^. As is shown in refer- 



4l| , the on-shell limit of equation (|2.5| ) recovers the result given in Q^'^^ is also 



known in the context of (9(a^) corrections to e^e annihilation into fermion pairs 



^4 separartes into contributions from the t- and u-channels and a t-u interference con- 
tribution (compare equation ( |E.12| )): 

e?r"(s;si2,s34) = Gl + g^ + GT 



Ga 



8 ■S12S34 
1 

8S12S34 



- + 2s(Si2 + S34) - 8S12S34 + 4si2S34(s-Si2-S34)£b 
D 



^ + 4s(Si2 + S34) + 16 S S12S34 (1 

3 V S-S12-S34, 



B 



[2.6) 
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If considered separately, the contributions Ql from the t-channel, Q2 from the u-channel, 
and from the t-u interference violate unitarity. As can be seen from equations (|2.4|) 
and ( ^.61 ), every contribution yields a total cross-section behavior proportional to s as 
s oo. Thus each single contribution violates the Froissart limit [0, and only sum- 



mation of the contributions Ql, Q^i Ga-, corresponding to the gauge invariant set of 



Feynman diagrams shown in figure 2.1, delivers a reasonable result with proper high en- 



ergy behavior. This is a standard result of Yang-Mills theories [Q. It is noteworthy that, 
adding the three contributions from t-channel, u-channel, and t-u-interference, the factors 
l/('Si2S34) contained in the individual kinematical functions given in equation ( p. 61) cancel. 
Consequently, the presence of the overall damping factor (si2S34)/s^ in dPa^ / {dsi2ds^4) 
should be considered a result of the gauge cancellations ensuring the unitarity of . 

A graphical representation of \/\/ s-Q^^'^ is given in plots a) and b) of figure |2]^ for two 
different ranges of the variables 1 — S12/ s and 1 — S34/S. As \fX/ s-Q*^^ is a dimensionless 
quantity, it only depends on S12/S and S34/S. From plot c) of figure ^]2| one can see how 
the double-resonance behavior entering the cross-section via the product pz(si2) Pz{.S2,i) 
of Breit-Wigner densities dominates the character of the twofold differential cross-sec- 
tion d?a^ / {dsi2ds'i^): Two intersecting resonances are seen in plot c) whereas the behav- 
ior of ^4'''" in dependence on S12 and S34 is rather suppressed. Final cross-section results 
are obtained from equation by twofold numerical integration^. 

To be able to compare the on-shell with the off-shell case, a "total" cross-section is intro- 
duced via 



a [s 



^tot(^) - 2-BR{l)-BR{2) ' ^^■'^^ 

which, in the narrow width limit, corresponds to the on-shell pair cross-section. The 
effect of the finite width on the total cross-section ( p.7| ) is illustrated in figure 2.3 . 
The off-shellness of the Z^ bosons causes the characteristic smearing of the excitation 
curve with a dumped peak and stretched tails. The change in the cross-section induced 
by the finite width is of the order T z/Mz and ranges from —7% at the cross-section peak 
to +5% at 1 TeV. The peak is shifted by approximately +4.5 GeV. Using a constant 
width, i.e. z{s) — > M^F^, corresponds to a redefinition of the Z^ boson mass ||60[| : 
Mz ^ Mz = Mz + \T\/Mz, resulting in Mz^Mz + 34 MeV. 

In the case of W pair production, gauge violating imaginary parts from the introduction 
of the finite width are only properly cancelled by imaginary parts from the trilinear vertex 



correction |62[. Such gauge invariance problems due to the naive introduction of the 
finite boson width into the propagators are not present in pair production. Possible 
gauge violations from the introduction of the Z° width as an s-dependent or s-independent 
quantity are phenomenologically irrelevant and may be trivially resolved by adjusting the 
definition of the width. The absence of the more serious gauge violations in off-shell 
Z^ pair production is understood from the fact that gauge cancellations already happen 
at the level of cxf because of the identical coupling structure of the t- and the u-channel 
contributions. 



^AU numerical programs for this thesis were written in Fortran. The author is grateful to D. Bardin 
for providing the integration routines SIMP and FDSIMP together with routines for the evaluation of 
Dilogarithm and Trilogarithm functions. 
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Figure 2.2: Plots a) & b): Behavior of the dimensionless kinematical function \/X/ s ■ Gl^"" 
for the two different ranges < 1—Sij/s < l—Sij^min/ s o,nd 0.15 < l—Sy/ s < 0.995. The var- 
iables I—S12/S and 1— S34/ s were chosen for the sake of clearer graphical presentation. Note 
that the used graphics package PAW / |5^/ has problems at the lower edges of phase space 
because of the very steep gradient of the plotted function. The flat region in the foreground 
represents zero level and lies outside the phase space. The edge of phase space is clearly 
seen from plot b). Plot c): Multiplication of VX/ s-Ql"^^ with s- pz{su)pz{s34) at = 
200 GeV emphasizes the dominance of the double-resonance character of process 
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Figure 2.3: The total pair production cross-section crj^^^{s) for process ( [ i . j| j and the 
effect of the finite width. The numerical input for the Z^ pair production figures is: 
6"/, = 1.16639 X 10-5 GeV-^, Mz = 91.173 GeV, TziMz) = 2.487 GeV, sin^^^ = 0.2325. 
Final state fermion masses are neglected. The numerical precision is estimated to be better 
than O.OlJoQ. 
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Chapter 3 

0{a) Initial State Radiation 



In e+e~ annihilation, initial state radiation (ISR) is known to represent the bulk of the 
radiative corrections. The 0{a) 'amputated' Feynman diagrams for initial state brems- 
strahlung to process ( p. .11) are shown in figure 3J., the corresponding virtual initial state 
correction diagrams are given in figure External leg self energies are absorbed into on- 
shell renormalization. For both the t- and the u-channel, bosons have to be attached 
to the two unconnected vertices of each 'amputated' Feynman diagram. 
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The computation of the virtual 0{a) cross-section starts from the 2 — > 4 particle phase 
space given in equation (|D.1|) , appendix p.l| . A cross-section differential in Si2 and S34 is 
obtained by fivefold analytical angular integration of the interference between the Born 
and the virtual matrix elements over the boson scattering angle and over the fermion decay 
azimuthal and polar angles in the corresponding boson rest frames. The computation of 
0{a) initial state bremsstrahlung starts from the 2 — >■ 5 particle phase space given in 
equation ( p.7| ), appendix p.2| . Then, a sevenfold analytical angular integration of the 
bremsstrahlung matrix element squared is carried out over the photon scattering angle, 
the boson scattering azimuthal and polar angles in the two boson rest frame, and over the 
fermion decay azimuthal and polar angles in the corresponding boson rest frames. After 
these five- and sevenfold angular integrations, the differential cross-section for the pair 
production process ( p..l| ) including 0{a) ISR with soft photon exponentiation is given by 



2^ISR 



dsi2 dsu 



(VsT2+V^«34)^ 



ds' r 1 

— pzisu) Pzisu)-2BR{l)BRi2)- P^^-'S + H (3.1) 

S L J 



with /3e = 2Q;/7r [ln(s/mg) — 1] and v = 1 — s'/s. This is exactly the result derived in ap- 
pendix |K^ and finally given in equation ([E.68|) , except that the sequence of integration has 
been changed to match the 2— i>5 phase space parametrization described in equation ( p.9| ), 
appendix p.2| . The soft -|- virtual and hard photonic parts S and H in equation ( |3.1| ) were 
calculated analytically. Both S and Ti. separate into a universal part with the Born cross- 
section at the reduced center of mass energy squared s' factorizing and a non-universal 
part: 



S{s,s';si2,su) 
H{s, s'; S12, S34) 



l + Si(s] 



•0-4 (S ; S12, 534) + 0-§[S ■,Si2,S34) 



ifi(s, s')-o-f (s; S12, S34) + o-^(s, s'; S12, S34) 



(3.2) 



Universal Part 



Non— universal Part 



The 0{a) radiators 5*1 and Hi were explicitly derived in the full calculation. They emerged 
from the virtual and bremsstrahlung matrix elements after five- and sevenfold angular 



integration. Their derivation is outlined in appendix E.5 



^1 
Hi 



a 

TT 



y 



1+^ 

s 



(3.3) 



(9(q;^) corrections may be taken into account in the universal corrections by adding the 
radiators S2 and H2 to Si and Hi respectively. 5*2 and H2 are e.g. found in reference pB| . 
The factorizing parts of the cross-sections S and Ti are called "universal" , because they 
are independent of the process topology, and especially of the initial state topology of 
the radiation proces^ In fact, the universal radiators Si and Hi are also found in the 
literature on fermion pair production in e^e~ annihilation, e.g. in reference It must 



be stressed that, although it was expected to find a separation of the cross-section into 



4n this statement it is implicitly assumed that the considered initial state topologies have unit charge 
initial state fermion lines connecting the initial state e"*" and e~. 
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universal and non-universal contributions, the working out of the separation was a techni- 
cally difficult task, especially for the t-u interference contribution. The virtual-|-soft and 
hard, non-universal, non-factorizable contributions in equation (|3.2D are given by 



as{g]Si2, S34j = - ■ g^^, ' {Lt + Re) " f^i,non«™ l^' ! ^12, S34J 
(Tnis, S'; S12, S34) = - ■ ^^^J^^^ te + ^e) ■ (^inonuniis, s' ] S12, S34) (3.4) 

with crYnonuni ^fnonuni derived in appendix ^. The full expressions are given in 
equation ( [E.20D , appendix |li].2| and in equation ( [ti].32D , appendix |E.3[ Equations ( |E.2CI| ) 



and ( [E.32|) represent a main result of this thesis, but will not be rewritten here because of 
their considerable length. The result ( ^.41 ) is new even in its limit for the case of on-shell 
bosons. It is seen from equations ( |E.2CI| ) and (|E.32| ) that the non-universal contri- 
butions contain many dilogarithm and trilogarithm functions and are rather involved. 
One should note that, in equation ( p.4| ), the virtual+soft non-universal contribution is 
evaluated at s' instead of s. The reason for this treatment is to avoid an unphysical, 
5 distribution-like concentration of non-universal virtual-|-soft ISR corrections at zero 
radiative energy loss. One of the reasons for the appearance of non-universal ISR con- 
tributions is the angular dependence of t- and u-channel propagators. For both virtual 
and bremsstrahlung corrections the angular behavior of t- and u-channel propagators is 
much more complicated than the angular behavior of s-channel propagators (see appen- 
dices p.2| , [E.2| , and |E.3| for details). Total cross-sections are computed from equation ( |3.1|) 



by threefold numerical integration over s', S34, and S12. 

Several other remarkable properties of equation ( ^.2| ) are worth mention. The leading ISR 
corrections to the cross-section with the mass singularities j3e factorize and are contained 
in the universal cross-section contributions. This is true not only for the overall cross- 
section, but also for the individual contributions from the t- and the u-channel and from 
their interference. As one would expect, non- universal contributions do not contain the 
mass singularity (3e- As a consequence, non-universal ISR corrections are suppressed with 
respect to the universal ones. The numerical influence of universal and complete ISR cor- 
rections on the cross-section for process ( |1.1| ) is presented in figure |37^ . Some discussion 



of the results presented in figure |3]^ is in order. Around threshold, in the energy range 
of LEP2, universal ISR peaks at large negative contributions with a relative correction of 
almost -30% at the pair threshold itself. Below A/i=280 GeV, ISR radiatively reduces 
the Z° pair cross-section. At higher energies ISR corrections are positive and develop a 
radiative tail increasing the Born off-shell cross-section by +14% at a/s=1 TeV. These 
patterns are familiar from the Z*^ pole, which is not surprising, because single Z° pro- 
duction is a resonant process too, and because the radiators Si and Hi describe the 
complete 0{q) ISR at the Z° pole |64|. The radiative tail, already present in universal 



ISR, is further pronounced by non-universal ISR, whereas non-universal corrections are 
hardly noticeable around threshold and therefore irrelevant in the LEP2 energy regime. 
The explanation for the radiative tail in the universal corrections is that radiative reduc- 
tion of s to s' increases the relevant Born cross-section erf (s', S12S34) (see equations (|2.4|) 
and ( 3.2|) ) in the region of negative slope of af. The additional corrections due to non- 
universal contributions range from 0.3% below threshold up to 2.5% at 1 TeV. In the 
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Figure 3.3: The total off-shell pair production cross-section (Jtot{s)- The solid line 
represents the Born cross-section crj^oti^) figure \2.$ The dash-dotted line gives the 

cross-section with universal ISR corrections only, i.e. and a^j are set to zero. The 
dotted line gives the cross-section with complete ISR from eq. i\3. In the inset, the 
relative deviations of the universally and completely ISR corrected cross-sections from the 
off-shell Born cross-section crj^oti^) ^''"^ given. The numerical input for the ISR corrections 
is: a = 1/137.0359895, me = 0.51099906 MeV . The numerical precision is around 0.1%q. 



energy range of LEP2, non-universal corrections do not exceed 0.4%. Comparison of fig- 
ures 2^ and 3^ shows that the high energy tail due to ISR is much more pronounced 
than the one due to the off-shellness of the bosons. It should also be realized that ISR 
shifts the cross-section maximum by 9 GeV from y^=211 GeV to i/i=220 GeV. 
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Vs [GeV] 


<t [pb] 


ISR.uni r i i 

(^tot [pb] 


(^tot [pb] 


130.0 


0.0015 


0.0013 


0.0013 


150.0 


0.0067 


0.0055 


0.0056 


165.0 


0.0215 


0.0175 


0.0176 


180.0 


0.1657 


0.1227 


0.1231 


190.0 


0.8971 


0.6995 


0.7021 


200.0 


1.1397 


0.9641 


0.9678 


212.0 


1.1898 


1.0629 


1.0671 


230.0 


1.1316 


1.0611 


1.0657 


260.0 


0.9799 


0.9613 


0.9663 


300.0 


0.8099 


0.8211 


0.8264 


400.0 


0.5503 


0.5811 


0.5866 


500.0 


0.4074 


0.4399 


0.4454 


600.0 


0.3171 


0.3481 


0.3534 


800.0 


0.2111 


0.2374 


0.2422 


1000.0 


0.1524 


0.1744 


0.1787 



Table 3.1: Numerical values for the off-shell pair production cross-section at tree level, 
with universal ISR, and with complete ISR as obtained from equations ( p^. i[ ) and i[ j by 
numerical integration. 



Investigating the unitarity for the ISR corrected cross-section, one observes from equa- 
tions ( |3.1| ) to ( |3.3| ) that the 1/s behavior observed in the Born cross-section is altered by 
factors ln(s/mg) originating from the mass singularity jSe- Therefore the unitary behavior 
due to the gauge cancellations discussed in chapter |^ is carried through to the universal 
part of the ISR corrections. For s — s>oo , the universally ISR corrected cross-section be- 
haves better than (ln^s)/s and thus approaches zero. As is seen from table O, the 



absolute value of the non-universal corrections reaches a maximum around y^=450 GeV 
and then decreases with increasing ^/s. This already hints at the unitary behavior of 
the non-universal ISR corrections which is ensured by the overall damping factor already 
discussed in chapter |]. However, an important difference between the appearance of 
the factor (si2S34)/s^ in the tree level and in the non-universally ISR corrected cross- 
section exists: All ISR cross-section contributions, namely the ISR t-channel, the ISR u- 
channel, and the ISR t-u interference contributions, contain this screening factor, whereas 
at tree level its appearance is only a result of the gauge cancellations. Thus the screening 
factor appears separately in all the following non-universal cross-section parts, 

/ \ / \ t,u,tu ^12 -534 /„ r-\ 

P{S12) p{S3i) ■ (T^Jj , (3.5) 

and it ensures the unitary behavior, i.e. the smallness of non-universal ISR corrections, 
at the level of individual contributions. This is more than one could expect, because 
cross-section unitarity only requires a unitarity-conserving mechanism for the sum of all 
contributions. 
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Figure 3.4: The generic two-boson production and decay Feynman diagram for e'^e anni- 
hilation. V stands for any vector boson, the fi, i = 1, 2, 3, 4 represent the decay fermions. 



To summarize this chapter, the previously unknown full initial state QED corrections to 
Z° pair production were calculated and separated into a universal leading logarithm and 
a non-universal part. It was illustrated how high energy unitarity is conserved for the 
radiatively corrected cross-section. Complementary results for W^W~ pair production 
were reported in |^ where, instead of a u-channel amplitude two s-channel amplitudes 
containing the W^W~ and W~^W~'y vertices are present. Taking these results into 
account, one can claim knowledge of all initial state QED corrections for the whole class 
of four-fermion production via two-boson production processes in e+e" annihilation. The 
corresponding generic two-boson production and decay Feynman diagram is shown in 
figure 
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Chapter 4 



Generalization to Neutral Current 
Pair Production 



Four-fermion production via massive boson pairs is easily accessible through experimental 
cuts on fermion pair invariant masses. But still one would like to know not only the 
four-fermion cross-section proceeding via pair production, but the full four-fermion 
generalization 

e+e- ^ /1/1/2/2 , /i^/2, U^e^M (4.1) 



of process ( |1 . 1| ) . Process ( [4.1| ) is described by 24 electroweak Feynman diagrams If 
/i and /2 are quarks, QCD contributions from 8 gluon exchange diagrams must be taken 
into account , and if massive fermions are produced, associated Higgs production must 
not be forgotten. In the electroweak case, if /i and /2 belong to the same isospin doublet, 
both charged current (CC) and neutral current (NC) processes contribute to reaction ([4.1D. 



The various aspects of the CC reactions have been studied by many authors 



68| , |69| , |fO| , |71| . Interferences between CC and NC contributions should be small due to 
different invariant masses favored by the production of fermion pairs via boson decays. 
Born level semi-analytical calculations for the NC case are available, and the inclusion of 
initial state radiation (ISR) in a leading logarithm approximation is unproblematic [^, 



65| . In figure ^]T|, all NC Feynman diagrams of the so-called Conversion type are depicted. 



As is seen from table |4.1| , the Conversion Feynman diagrams strongly dominate the cross- 




Figure 4.1: The Born level Feynman diagrams for neutral current pair production. Left: 
t-channel. Right: u-channel. 
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v/i [GeV] 


^Nc [fb] 

All TV T /"I 1 • 

All Mu diagrams 


U^U bb [n 1 

(^conv [fb] 
NC Conversion diagrams 


91.0 


613.832 


278.080 


141.0 


94.469 


94.477 


191.0 


53.252 


53.231 


241.0 


38.610 


38.613 


291.0 


28.396 


28.422 


341.0 


22.128 


22.161 


391.0 


17.930 


17.970 


441.0 


14.946 


14.982 


491.0 


12.713 


12.742 


591.0 


9.620 


9.570 


691.0 


7.551 


7.603 


791.0 


5.935 


6.187 


891.0 


5.166 


5.165 


991.0 


4.237 


4.382 



Table 4.1: Neutral current (NC) fi^fi^bb production cross-sections. By 'NC" it is empha- 
sized that only amplitudes containing intermediate Z° bosons or photons are accounted 
for. Feynman diagrams containing the Higgs boson are excluded. No cuts are applied. 
The author is grateful to A. Leike for supplying the Fortran code 4fAN ^6d{] for the 
production of the results in this table. 



section at energies above the Z° mass. Interferences between Conversion type and other 
NC diagrams are small and negative. Therefore it is interesting to investigate the process 

e+e- ^ (Z°Z0,Z%,77) - hhf2Ul) , /it^/s , f^^e^X (4.2) 



with full ISR, proceeding via the tree level diagrams shown in figure |4.1j . The exchange 
photons to be taken into account have no resonance behavior. Consequently, contrary to 
the double-resonance case treated in chapter ^, it is not advantageous to express initial 
state couplings and boson propagators in terms of Breit-Wigner densities as given in 
equation ( p.2| ). Since process ( [4.2| ) can be described by several sets of t-channel, u-channel 
and t-u interference contributions, it is appropriate to replace equation (|2.1|) by 

with the twofold differential cross-section factorizing into a phase space factor, a function 
C of couplings and propagators, and the kinematical function Ql~^^ already known from 
equation (|2.5|), 



dsi2 dsu 71" 



•C(/i,si2;/2,S34)6^r"(s;si2,S34) • (4.4) 
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As ^4^" is, in the language of Feynman diagrams, related to a pair of t- and u-channel 
diagrams, the function C is obtained by summing over all 16 such pairs that can be formed 
from the eight diagrams of figure [4.1| . This results in 

2 Si2 Su 



C(/l,Si2;/2,S34) 



E 



(67r2)2 

1 



B„Bj,Bk,Bi='r,Z 



DbXsi2) Db,{s^,) D%^{S^2) D%{s^,) 



■X 



Le+e-Bi Le+e-Bk Le+e-Bj Le+e-Bi + Re+e- Bi Re+e- B,, Re+e-Bj Re+e- Bi 
^fJlB,LfJ,Bk + RfihB.RhhBk -Neifl) Np{Bi,Bk,mi,Si2) X 



^M2BjLfJ^Bi + R 



hhBj RhhBi 



■Ne{f2)Np{B„Bi,m2,Su] 



X 



(4.5) 



where left- and right-handed couplings are used as defined in equations (p.27| ) and 
The propagator denominators are given by 



DBiSi 



Mb + iy^TB{s 



(4.6) 



Sij Tz/Mz 



with their complex conjugates D*^. = = for the photon and Tz{sij) - 
which is a good approximation around the peak and above. The color factor Nc{f) is 
unity for leptons and three for quarks. The phase space factor 



Np{Bi, Bk, m, s) 



for Bi^ 'J or Bk ^ j 



1 - Am'^/s ■ (1 + 2mVs) for Bi = Bk 



(4.7) 
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correctly takes into account that fermion masses may be neglected if the corresponding 
fermion pair couples to a Z° resonance but may not if it only couples to a photon. In 
the latter case, finite cross-section contributions arise from close to the lower boundary of 
the invariant masses Sij. It is readily verified that the above generalization of the pair 
results from chapter ^ is formally carried out through the replacement 



pzisu) Pz{su)-2BR{1)BR{2) ■ 



{G,Ml 



{Lt + Rt 



C{fi,s,2-J2,su) . (4.8) 



-2 



The input for numerical calculations was chosen as = 1.16639-10 ^ GeV 
91.173 GeV, Tz = 2.4971 GeV, Mw = 80.22 GeV, and = 0.51099906 MeV. For 
the electroweak mixing angle, the effective value sin^^^ = 0.2325 was chosen. The input 
for the left- and right-handed couplings was Qe = — I, Tw\ ~ ~\i ^^"^ ^ ~ V^vra with 
a = \/2My^ sm^Ow G ^/ Ti . For the QED corrections, because they are mainly due to soft 
photons, the Thomson limit a = 1/137.0359895 [10| was used. These parameter choices 



enabled straightforward comparisons with the results in reference ||65|. For several center 
of mass energies in the range of LEP2 and above and for different final state fermion pair 
combinations, tree level results for the full set of Conversion diagrams were compared 
to the appropriate contribution as obtained from the Fortran code 4f AN used in ref- 
65 1 . Agreement was found to be at the level of the numerical precision of the 



erence 



integration routine. 
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Taking ISR into account, means to implement the replacement ( [4 .81 ) into formula ( |3.1| ), 
resulting in 



f42 ISR 
" '^Conv 

ds 12 ds^i 



ds' 



C(/i, S12; /2, S34) ■ iSeV'^" Sconv + 'Hcon 

S L 



(4.9) 



for the twofold differential 0{a) ISR corrected cross-section with soft photon exponen- 
tiation for the NC pair production process ( |4.2| ). Recalling that the factor (G^Af|) /8- 
(Lg + _Rg) is now absorbed in C and defining 



TT S 



2 ^4^"('5; ■S12, S34) 



^Coral-^; ■512, -534) 



one finds from equation ( p.4[ ) and in analogy to equations ( |3.2|) and ( p^ 



(4.10) 



Sconv{s, s'; S12, S34) 



l + 5i(s) ■Or^o„^(s';Si2,S34) + (^5,Con^ (^'^ ^12. ^34) 



Conn 



S,s';Si2,S34) = Hi{s,s')-aconv{s''^Sl2,S34) + "^H Co™(s, s'; ^12, S34) 



Universal Part 



Non— universal Part 



(4.11) 



with the radiators Si and Hi given in equation ( p.3| ) and 



f^i?,Con,;lS' « ; S34, 



« '^rnonum('5'; ■S12, S34) 



TT 



pi? 



TT 



,S, s'] S12, S34) 



(4.12) 



(^Inonuni and are found in equations ( |E.20|) and (|E.32| ) respectively. 



As it yields a final state of experimental importance for Higgs searches and as it is also 
representative for semi-leptonic final states, the following presentation of results will focus 
on the reaction 



e e 



(4.13) 



which has already received considerable attention in the literature |75|, 65|. In fig- 
ure cross-sections for reaction ( |4.13| ) as obtained from equations ( [4. 3D and ( |4.9| ) are 
presented. The cross-section correction due to universal ISR is approximately 12% at 
yi=130 GeV, increases to 18% below the Z° pair threshold, decreases to 13% in the 
region of the Z° pair- induced cross-section "bump" and then steadily increases to 21% at 
600 GeV. The additional relative correction from non-universal ISR corrections increases 
from 9 %o at 130 GeV to 4.2% at 600 GeV. These numbers can be extracted from table HT2 



where numerical values for the cross-sections presented in figure |4.2| are given. The nu- 
merical precision of the values in table is around 0.1 ^/ooQ With respect to pair 
production the relative corrections of both universal and non-universal ISR are enhanced. 

^In the following, numerical precisions will not be quoted any more. They are around 0.1 %o or better 
throughout this thesis. 
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Figure 4.2: The NC pair production cross-section for process d^.lSj ). The solid line repre- 
sents the Born cross-section, the dash-dotted line includes universal, and the dotted line 
includes all ISR corrections. In the inset the full NC pair cross-section is compared to the 
contributions from and photon pair production. The muon and b quark masses were 
toA;en as m^ = 0.105658389 GeV andmb = 4:.3 GeV ^ 



Figure 4.3: The effect of cuts of 2-T z and h-T z around the Z mass Mz on the NC pair 
('All Graphs') and Z^ pair ('ZZ Graphs') cross-sections. Tz = 2.4971 GeV. The cuts 
were applied to both fermion pair invariant masses. All curves show cross-sections for 
reaction ( U.lSj ) with universal ISR corrections. 
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[GeV] 


%\-.5 [fb] 


ISR.uni rn 1 


Vjr^'- [fb] 


ISR.uni rn 1 

[fb] 

„/ / „ \ n o 
s 1 s dL u.y 


vjr^'- [fb] 

s / s ^ u.y 


130.0 


131.4968 


147.4554 


148.7304 


111.5411 


112.2462 


150.0 


84.3863 


97.9860 


99.2177 


71.1558 


71.7510 


165.0 


65.8298 


77.5651 


78.7483 


55.3251 


55.8606 


180.0 


55.0144 


64.8390 


65.9749 


45.9486 


46.4371 


190.0 


56.4100 


63.8071 


64.9309 


46.6443 


47.1260 


200.0 


54.0441 


60.9149 


62.0168 


44.9442 


45.4106 


212.0 


49.8441 


56.4425 


57.5120 


41.5066 


41.9482 


230.0 


43.7651 


49.9910 


51.0124 


36.4290 


36.8391 


260.0 


35.7550 


41.3205 


42.2653 


29.7064 


30.0696 


300.0 


28.3717 


33.1358 


33.9928 


23.5131 


23.8300 


400.0 


18.1597 


21.5427 




14.9760 




500.0 


12.9578 


15.5310 


16.1086 


10.6475 


10.8406 


600.0 


9.8428 


11.8958 


12.3895 


8.0646 


8.2251 


800.0 


6.3587 


7.7901 




5.1868 




1000.0 


4.5111 


5.5883 




3.6671 





Table 4.2: Numerical values for the cross-section of process ([j.lSJ at tree level, with 
universal ISR, and with complete ISR as obtained from equations ( U-!^ to The 
last two columns give ISR corrected cross-section values after cutting on s' as indicated. 



Contrary to Z° pair production, ISR corrections are always positive which is easily ex- 
plained by the negative slope of the cross-section curve for all center of mass energies a/s- 
Thus a radiative reduction of s to s' enhances the cross-section in a way very similar to 
the radiative tail in pair production. It must be stressed that the absolute value of 
the non-universal corrections to reaction ( [4.13| ) steadily decreases with energy, which is 
expected as an effect of the screening property. It is seen from figure [4.2| that, as expected, 
ISR smears out the pair threshold behavior. In the inset of figure the full NC pair 
production is compared to the contributions due to Z° and photon pair production for 
the universally ISR corrected case. Clearly, the main part of the cross-section is due to 
contributions which contain both bosons and photons. 

Figure ^]3| shows how the Z^ pair production cross-section is restored from the full set 
of Conversion diagrams through invariant mass cuts on Si2 and S34. Also seen is the 
influence of the applied cuts on the Z^ pair cross-section. 

In figure [4.4| it is shown how the cross-section for reaction ( [4.13|) evolves when invariant 
mass cuts on the /i"*"//" and the bb pairs are tightened. Already cuts of 50 GeV on both 
invariant pair masses yield an s-dependence of the cross-section very much resembling 
the double- resonant excitation curve for Z° pair production. The cuts |M^+^- — Mz\ < 
5Tz and M^^ > 60 GeV represented by the dash-dotted line in the figure are of interest for 
Higgs searches, as the current Higgs mass limit lies around 60 GeV . In these searches. 
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Figure 4.4: Evolution of the cross-section for reaction ( 4-1!^ for different cuts on the 
H~ and bb invariant pair masses. All curves show cross-sections with universal ISR 
corrections. 



associated Higgs production is used with the bb pair couphng to the Higgs boson. 



From figure ^]5| one can see how cross-sections behave if cuts are apphed to the reduced 
center of mass energy squared s', i. e. to the invariant mass squared of the four-fermion 
system. In experimental analyses such cuts are very commonly used to reduce back- 



grounds. Numerical values after cuts on s' are given in table |4.2|. Both plots of figure ^75 



show cross-sections for the final state /i^ bbi^). The upper plot is for the pair reac- 
tion ( |1 . ![ ) , the lower one for process ( |4.13| ). The main windows present the development 



of universally ISR corrected cross-sections when the s' cut is tightened. The insets, where 
the ratio of the cross-section with complete over the cross-section with universal ISR 
is plotted, give the additional relative cross-section contribution from non-universal ISR 
correction. From the upper plot for e+e~ — > (Z^Z^) —>■ fi^fi^bb one can make three observa- 
tions. First of all, one notices that the cross-section with ISR is dominated by rather soft 
radiation. Secondly, hard photon ISRQ is strongly suppressed below the double-resonance 
peak. Only above the peak, such hard photon events represent a significant cross-sec- 

^^By hard photon radiation the author denotes the radiation of photons strongly reducing s. This 
"strong reduction" may be several per cent or more than 10%, depending on the context. 
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Figure 4.5: Effect of s' cuts on the fj,~^fi' bb {-y) cross-section for final state production via 
process ^Tlj ) (upper plot) and process ( U-^^ (lower plot). The main windows present uni- 
versally ISR corrected cross- sections, the insets show ratios of completely over universally 
corrected cross-sections. No cuts are applied to Si2 and S34. 
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tion contribution. Both these results are analogous to similar ones at the single Z° peak 
at LEPl [0. This is, because both processes, single and double production, are 



resonant processes with Breit-Wigner factors 
\s-Ml: + t^Tz{s)\^ 

limiting hard radiation in the peak region. In addition, both calculations use the same 
0{a) radiators Si and Hi. The third observation, made from the inset of the upper plot, 
is that non-universal ISR tends to be harder than universal ISR so that the ratio of the 
cross-section corrected with complete ISR over the universally corrected cross-section de- 
creases with a cut on s'. In other words, the non-universal ISR corrections in per cent 
decrease with a cut on s'. This may be understood as follows: 0{a) non-universal correc- 
tions are infrared finite, whereas the infrared divergent universal corrections underwent 
a resummation, namely the soft photon exponentiation. Thus the universal corrections 
contain important soft resummed parts, and the non- universal corrections don't. The 
opening angle between the two curves in the inset reflects the increasing importance of 
hard radiation at higher energies. Around the double peak and below however, the 
ratio of the completely over the universally corrected cross-section is unaffected by the 
s' cut. This means that, below ^/s ~ 200 GeV, no considerable excess of hard photon 
radiation due to non-universal corrections exists. As for universal ISR, where deviations 
of the cross-section with the cut s'/s > 0.9 from the cross-section without cut are only 
seen above -\/i~200 GeV too, harder radiation is disfavored by the steep rise of the tree 
level cross-section around threshold. 



Turning to the lower plot of figure which is for reaction ( |4.13| ), patterns similar to the 



upper plot of the figure are encountered. There is, however, no disincentive for hard radia- 
tion due to a steep rise of the cross- sect ion. Therefore, the ISR corrected cross-sections in 
the main window are lowered throughout the presented energy region by all cuts on s'. 
The approach between the curves for no cut and for s'/s > 0.5 around 150 GeV must be 
attributed to the inhibition of a reduction of s' below the threshold for single-resonant 



Z°7 — >■ ^"hh production. From the inset of the lower plot of figure f4.5| one may draw 
the same conclusions as for the inset of the upper plot, except that there is no disturbance 
of hard radiation due to the rise of the cross-section in the double Z^ resonance region. 

To conclude the discussion of figure |4.5| , it is mentioned that the integration of the in- 
volved non-universal cross-section contributions is very CPU time consuming. This means 
that, contrary to what would be desirable, it is not possible to mass-produce results with 
non-universal ISR and all kinds of cuts. Still, a more thorough analysis including an 
investigation of the spectrum of radiated photons is under way |]72|. 



Figure |4.6| shows cross-sections for process (|4.2| ) with different final states. The upper 
plot of the figure shows final states containing hb pairs, the lower plot the smaller cross- 
sections for final states with /i"*7i~ pairs. Cuts were applied to require y^si2>5 GeV and 
a/ -5 34 > 20 GeV. Comparing the final states uuhh and ddhh, it is interesting to observe how 
the different u and d quark electric charges affect the cross-section below threshold: There, 
the ddhh cross-section is much smaller than the uuhh cross-section . Similarly interesting 
is the comparison between fi~^fi~hb and v^v^hb. The v^v^hh cross-section has a shape very 
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Figure 4.6: Universally ISR corrected cross-sections from process ( U-W for different four- 
fermion final states containing bb pairs (upper plot) or fi~^fi~ pairs (lower plot). 
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Fermion 






r 


u 


d 


c 


h 


Mass [GeV] 





0.105658389 


1.7771 


0.005 


0.01 


1.3 


4.3 



Table 4.3: T/ie fermion masses for the numerical computations presented in the figures. 



much like the Z*^ pair excitation curve due to the cuts and the weak coupling of b quarks 
to photons. The fi~^fi~bb cross-section on the other hand receives significant contributions 
below the pair threshold due to the muon's electric charge. Also rather interesting, but 
straightforwardly explained by the cuts and the photonic and weak couplings, are com- 
parisons between the final states fi~^fi~bb versus fi~^fi~cc and versus fi~^fi~ t~^t~ . 



The fermion masses for all numerical computations were taken from reference |]T0[ and 
are listed in table 0|. 



To conclude the chapter, figure \i.7\ presents the infiuence of final state fermion masses 
as taken into account in the semi-analytical approach. It is recalled that final state 
fermion masses are only taken into account through the factors Np{Bi, Bk,m, s). Fig- 
ure |4.7| shows, in the main window, the two cross-sections of reaction ( [4.2| ) into the final 
states fi~^fi~dd and fi~^fi~bb which do not differ in couplings, but only in final state quark 
masses (see table The J cross-section is larger because of its larger phase space 

due to nid < rrib. If now the invariant dd mass is cut at exactly two times the b quark 
mass, the phase spaces for the two final states have equal volume. However, a difference 
remains due to the factors Np for the quarks. While for the final state fi~^fi~dd 

iVp(7,7,mrf,S34) = y'l - 4m^/s34-(l + 2ml/ S34) ~ 1 



throughout the phase space cut at ^/su = 2 nib, because of the, compared to the cut, 
small d quark mass, this is different for the final state ^'^ix~bb. iVp(7, 7, rrife, S34) is zero at 
the phase space boundary ^/s^ = 2 and approaches its limiting value 1 to a precision 
of l%o only around 1/534 = 38 GeV. Thus the cross-section for the final state fi^fi^bb is 
slightly smaller than the one for fi^fi^dd cut at ^/s34 = 2mb. This is seen from the inset 
of figure where the ratio of cr^+^~^^ with the cut ^/su>2mb over a^+^-bi is given. 
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Figure 4.7: The effects of masses in the cross-sections for (i^jji'dd and n~^n~bb production 
according to reaction l \4-3i )- The main window shows the cross-section enhancement due 
to the larger phase space for the production of lighter fermion pairs. The inset shows the 
small effect of the factors Npipf, 7, m, s), if phase spaces are made equally large by applying 
a cut of ^/s^ = 2mb to the invariant dd pair mass. A small bump is seen around the 
pair threshold. All curves are with universal ISR. 
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Chapter 5 

Summary, Conclusions, and Outlook 



Even around the pole, four-fermion production is an interesting topic, relevant not 
only for Higgs boson searches |Q but also for higher order corrections to the line 
shape Above the production thresholds for resonant gauge boson pair pro- 



duction, e~^e~ annihilation into four fermions is one of the most significant Standard Model 
processes. Thus, double-resonance production plays a very important role at LEP2 and 
other future e^e~ colliders, and the study of four-fermion final states will provide further 
understanding of gauge boson properties and may be the key to the mechanism of elec- 
troweak symmetry breaking via investigations of associated Higgs boson production. 

In the preceding chapters, the first complete calculation of initial state QED corrections 
(ISR) for neutral current double-resonance production in e"'"e~ annihilation was presented. 
A semi-analytical method was used to compute cross-sections and invariant mass distri- 
butions for and neutral current (NC) pair production processes (|1 . 1| ) and ( |4.2|) p]. These 
processes are interesting by themselves and also as background to W^W~ physics and 
Higgs searches. All angular degrees of freedom were integrated analytically. After the 
angular integrations, only two fermion pair invariant masses were left for numerical in- 
tegration with a straightforward self-adaptive Simpson algorithm. In the case of ISR, 
also the four-fermion invariant mass had to be integrated numerically. It was shown that 
boson off-shellness and ISR both yield considerable corrections to the cross-section, each 
of C(10%) and up to C(20%) for neutral current pair production. ISR corrections sepa- 
rate into universal parts with the Born cross-section factorizing and non-universal parts. 
Contrary to non-universal contributions, universal ISR parts contain the large logarithm 
ln(s/mg) as a result of collinear photon radiation from the initial state electrons. Thus 
universal ISR is enhanced with respect to non-universal ISR. Invariant mass cuts can be 
easily implemented to isolate the Z° pair signal from the full NC process. Such cuts do 
not only suppress NC pair production, but also irreducible backgrounds of the annihila- 
tion type. Two examples of annihilation type background FejTiman diagrams are given 
in figure |5]^ (compare reference p5|). 

The systematic errors of the presented calculation are physics dominated, because of the 
numerical integration's excellent precision. Errors arise via 



'''A similar calculation for W^W~ pair production had to face ambiguities of the definition of initial 
state radiation due to the charge transfer from the initial to the final state in the t-channel W^W~ pro- 
duction amplitude p3. 
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Figure 5.1: Two examples of annihilation type Feynman diagrams for singly resonant 



electroweak background to reaction ( l^.Sj ). 



The use of the ultrarelativistic approximation (URA) with neglect of terms of 
0{ml/s) and 0{m1/ Sij). Errors due to the URA are very small. This is evidently 
true for pair production. As all used fermion masses are at least one order of 
magnitude larger than the electron mass, and as lower invariant mass cuts of several 
GeV are quite commonly applied, the URA is a valid approximation also for the 
NC pair production process (|4.2| ). 

Neglect of final state fermion masses in the matrix elements. For light fermions, the 
effect on the total cross-section is small. In the NC pair case with exchange photons, 
all important mass effects are properly taken into account by the factors Np given in 



equation ([4.7| ). Mass effects are removed by even moderate invariant mass cuts. For 
some cases, such cuts are also desirable with respect to non-perturbative threshold 
effects. 

Neglect of final state QED corrections. Compared to ISR, final state QED (FSR) is 
known to be suppressed, but may become important at very high energies beyond 
the TeV level. Interference effects between ISR and FSR were found to be suppressed 
by O^Tb/Mb) in resonant pair production of bosons B WM. 



Neglect of true electroweak corrections. For on-shell pair production true elec- 
troweak effects were found to yield corrections around 1% to 2% This estimate 
may have to be increased to accommodate off-shell boson pair production, but no 
full calculation of electroweak one loop corrections for process (^TTj) is available at 
present. In reference [^, a strong dependence of electroweak corrections on the 
Z^ scattering angle was found for center of mass energies above 500 GeV. 

Theoretical and experimental uncertainties of the parameters entered into the nu- 
merical calculation. 
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In the following, it will be shown how the investigations presented in this thesis are em- 
bedded in the environment of four-fermion final state physics and in the framework of 
semi-analytical calculations dealing with four-fermion production in charged current and 
neutral current processes El], |50, H3I |75|, |6^, 0]. The semi-analytical method yields 



compact formulae for total and differential cross-sections at tree level and with leading ini- 
tial state QED corrections. At present, semi-analytical cross-sections are available which 
are differential in the two boson invariant mass and the boson scattering angle . These 
are easily numerically integrated to yield singly differential or total cross-sections. In this 
thesis, the semi-analytical method was used to present non-leading initial state QED cor- 



rections |iT, 1^ 



The screening property was found in these non-leading corrections. A 
complete initial state QED calculation for a double-resonant four-fermion production is 
documented in detail for the first time. 

In the presented calculation, some non-resonant annihilation type background contribu- 
tions, see figure |5]3, to process (^31 ), 



e e 



/1/1/2/2 



are not taken into account. But already for moderate cuts, the results given in this thesis 
represent the dominant contribution to process ([4.1|) . Process ( [4.1D was completely and 
semi-analytically calculated in references |E0|, ITSl, B5[, but only at tree level. However, as 



background contributions are small at LEP2 energies and above and as the leading uni- 
versal ISR corrections may be straightforwardly implemented by multiplication of the tree 
level cross-section with the radiators 5*1 and Hi (see equation (|37^ )), sufficiently precise 
results are available for process ([4.11). This is still true if associated Higgs production is 
considered, because interferences of the Higgs Feynman diagram are negligible and because 
and Hi describe the Higgs diagram's complete ISR correct ionsQ. Thus, for the leptonicQ 
and the important semi-leptonic channels of process ([4.1|) , the description obtained from 
this thesis and references 



Jo], |75|, [6^ is sufficient in view of future e~^e~ collider data. The 
treatment of four-quark final states requires the consideration of eight additional Feyn- 
man diagrams of the annihilation type with decaying gluons radiated from the final state 
quarks. Two examples of such diagrams are given in figure p.2| . These diagrams were 
taken into account in semi-analytical and in Monte Carlo calculations |76|. The 



annihilation type QCD background diagrams are described by the same semi-analytical 
kinematical functions as the electroweak annihilation type background diagrams. The dif- 
ferences lie in coupling constants and color factors. Going a little further one realizes that 
hadron jets stemming from light quarks or gluons are experimentally not really distin- 
guishable so that one must incoherently add several contributions to obtain a prediction 
for an experimental signature. This means that new contributions with the virtual gluon 
from figure f).2\ decaying into two gluons instead of two quarks must be taken into account. 



^Generally, for Feynman diagrams of t he A nnihilation type with only one vertex in the initial state 
electron current (like the diagrams in figure 5A) the radiators Si and Hi d escr ibe the complete 0{a) ISR. 
As soon as Conversion type Feynman diagrams (like the diagrams in figure O) are involved, non-universal 
ISR contributions are present. 

■''It should be noted that final states with ^- and r-neutrinos can be calculated semi-analytically, 
but are not experimentally observable individually, because the signature from final states with electron 
neutrinos is identical. 
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Figure 5.2: Two examples of annihilation type Feynman diagrams for singly resonant 
QCD background to reaction U-W - The curled lines represent gluons. 



Recently, the reaction e'^e^ bb + 2 Jets was evaluated ||7g]. Another problem of the 
four-quark final states is related to so-called color reconnection effects which are mainly 
of interest for W^W~ pair production and the W mass measurement ||77|| . 

It is, at this point, interesting to compare the results for pair production with those 
of W^W~ pair production. Both calculations are fundamentally different in two regards. 
Firstly, the contributing amplitudes differ: t- and u-channel amplitudes for Z° pairs, but 
t- and s-channel amplitudes for W^W~ pairs (see e.g. reference |^^). Secondly, trilinear 
boson couplings contribute to Standard Model W^W~ pair production, but not to Stan- 
dard Model pair production. On the other hand there are, of course, patterns that are 
common to W^W~ and Z° pair production, such as the double-resonance behavior and 
the screening property. For four-fermion final states containing fermion-antifermion pairs 
from the same weak isospin doublet, double-resonant contributions arise from W^W~ as 
well as Z^ pair production, not to mention the many background contributions to both. 
The individual charged current and neutral current contributions to such final states have 
already been semi-analytically calculated. The interferences between CC and NC contri- 
butions, though small, require special treatment. This is due to different particle pairings 
from the decays of charged and neutral bosons which invalidates the unique assignment 
of the invariant masses S12 and S34 to fermion pairs. This problem is illustrated in equa- 



tion ( |5.1| ) for the double-resonant signals yielding a UUDD final state, where U and D 
denote generic particles with third weak isospin component +| and — | respectively: 



CC : 



e e 



I — > UD (inv. mass S12) 
U UD 



NC : 



e e 



r 

U DD 



[imr. mass S34J 
UU (inv. mass S12) 



(5.1) 



Imv. mass S34I 
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10 
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10 
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Table 5.1: Number of tree level Feynman diagrams contributing to the production of four- 
fermion final states via neutral bosons. 



The complete process e"'"e~ — > UUDD including the small CC-NC interferences is, at 
present, only calculated by Monte Carlo programs (see e.g. [0]). Similarly four-fermion 
final states containing electrons or electron neutrinos in the final state have not yet been 



calculated semi-analytically, but constitute no problem to Monte Carlo programs [|78| , [79 



The additional complexity of final states with electrons or electron neutrinos is due to 
their additional amplitudes with t-channel boson propagators. In table |5.1| , the number of 
Feynman diagrams contributing to four-fermion production via neutral boson exchange, 
i.e. the final states this thesis is concerned with, are listed. The simplest case (numbers in 
boldface) does not contain identical particles, electrons, electron neutrinos or particles 
from the same isospin doublet. Those final states were calculated semi-analytically and 
the calculations presented in this thesis describe the corresponding dominant contribu- 
tions for energies above the pole. For pairs of identical fermions, antisymmetrization 
of the final state as needed to satisfy the Pauli principle increases the number of Feynman 
diagrams (given in typewriter). Numbers in roman are for final states with electrons or 
electron neutrinos, and numbers in italics correspond to final states with particles from the 
same isospin doublet. Numerous Monte Carlo approaches [^, ^ ^ |73|, |7^, |7^, ^ 
treat four-fermion production in e"'"e~ annihilation. The scopes of the two independent 
approaches represented by the semi-analytical and the Monte Carlo method are different. 
The semi-analytical method provides formulae to strongly promote the understanding of 
four-fermion production processes. Monte Carlo programs on the other hand aim at com- 
plete sets of Feynman diagrams and are needed by experimentalists to produce input for 
their detector simulations and to apply all experimental cuts. Semi-analytical programs 
can only be run with quasi-experimental cuts. Unfortunately, Monte Carlo programs are 
often "black boxes" to the experimentalist. It becomes clear that semi-analytical and 
Monte Carlo calculations are nicely complementary. Keeping in mind that four-fermion 
final states are much richer in structure, much more interesting and much more compli- 
cated than fermion pair final states, the ultimate goal of semi-analytical calculations is 
rather obvious, namely to describe four-fermionic complexity by as much simplicity as 
possible. The semi-analytical method has been very successful in this respect, because a 
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e+e fi^fi bb 


e+e —>■ udfi 


e+e- ^ (^0^0,^%, 77) /i/i/s/sIt) 




25 


10 


8 (+56) 




325 


55 


36 (+532) 




4 


8 


1 (+2) 


Refs. 


[75, 65] 


|42, 81| 


This thesis 



Table 5.2: Examples for the number Ndiag of Feynman diagrams contributing to the pro- 
duction of four-fermion final states in e~^e~ annihilation. Also given are the numbers of 
interferences between Feynman diagrams (Ninter) cind the numbers of kinematical func- 
tions Ng needed to semi-analytically describe the process. The numbers in parentheses in 
the last column give the additional diagrams, interferences, and kinematical functions due 
to non-universal ISR as described in this thesis. 



description of many Feynman diagrams and their interferences by few kinematical func- 
tions (such as ^4"*"", see equation (|2.4| )) was achieved for several interesting four-fermion 
final states. Table |5.2| gives some examples. 

It is of fundamental theoretical interest to continue the systematic semi-analytical treat- 
ment of e~^e~ annihilation into four fermions. For the future it is planned to extend the 
work of this thesis to include a non-standard, heavy new neutral boson Z' and to evaluate 
the corresponding discovery limits at future e~^e~ colliders. Further a full calculation for 
process ( [4.1| ) with all ISR corrections and semi-analytical calculations of four-fermion final 
states with electrons and electron neutrinos are very desirable. It could also be interest- 
ing to extend the presented results to accommodate anomalous couplings between three 
neutral gauge bosons. Such couplings are not present in the Standard Model. Their in- 
clusion would require to take into account an additional Feynman diagram with s-channel 
e~^e~ annihilation into one neutral gauge boson which then decays into two neutral bosons 
that finally decay into the four final state fermions]^. It is worth mentioning that all ini- 
tial state QED corrections for such an inclusion are known from references [|T, |2| and 
from this thesis. Anomalous trilinear neutral gauge bosons couplings are discussed in 



reference [Q, but no deviations from the Standard Model were experimentally found so 
far |8^. Finally, it should be realized that the results presented in this thesis can be 
validated for proton-proton collisions at the Large Hadron Collider (LHC) by folding in 
parton distributions, changing coupling constants and introducing color factors. 



^Depending on cuts, one may also have to account for annihilation type background diagrams. 
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Appendix A 

Metric and Conventions 



Every calculation in four-dimensional space-time requires the choice of a metric which is 
uniquely defined by a metric tensor g^^. For this thesis, 
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V 



/i, = 0, 1, 2, 3 



(A.l) 



is chosen as metric tensor. This is just the negative of the Bjorken-Drell metric tensor |90 
A contravariant 4-vector is defined by 



(A.2) 



with time-like component w° and space-like components v. The corresponding covariant 
vector is given by 



gtMuV 



(A.3) 



where the convention to sum over repeated indices is adopted. The summation convention 
is used throughout this thesis unless stated otherwise. The scalar product of two 4- vectors 
is defined as the contraction with the metric tensor: 



v-w = Vfj^w 



V^Wa 



(A.4) 



It should be noticed that, with the above convention, time-like four- vectors have negative 
and space-like four-vectors have positive metric. The choice of metric also entails that 
the invariant mass squared of a system /C with 4-vector k is given by 



Ml = {ky 



The Dirac algebra is defined by 

75 



-k-k . 



(A.5) 



- i 7o7i7273 



(A.6) 



Finally, it should be mentioned that natural units c = h = 1 are used throughout this 
thesis. 
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Appendix B 

The Standard Model of Electroweak 
Interactions 



B.l The Standard Model Lagrangian 

Phenomenologically, the Standard Model is based on the empirical fermion content of 
nature given in table |B.1|, on the existence of the four vector bosons 7, W~^, W~ , and 



that mediate the electroweak force, and on the Gell-Mann-Nishijima relation ||84 

Qf 



iB.l] 



These empirical facts can be embedded into a theoretical description which is characterized 
by a Lagrangian that is invariant under SU{2)l x f/(l)y^ gauge transformations, 



SM 



Cg + Ch + Cf + C 



Yukawa 



(B.2) 



The Gauge Field Lagrangian 

The pure gauge field part Cq is given in terms of the massless gauge fields PV^ (a = 
1,2,3) and B^. The isotriplet M/^ corresponds to the three generators 7^ of the isospin 
group SU{2)l, belongs to the adjoint representation of SU{2), and couples to the left- 
handed doublets only. corresponds to the generator —Yw/2 of the U{1) subgroup 
of gauge transformations and couples to hypercharge. Adopting the conventions from 
appendix ^ one finds 

^G = -\w;;,W'''^''' -^B.^B^^'' . (B.3) 
Summation is also assumed over a, and the field strength tensors are given by 

B^"" = d^B'^-d'^B'' , (B.4) 

where eabc are the SU{2) structure constants and (71,(72 denote the U{1) and SU{2) cou- 
pling constants. The non-abelian structure of the gauge group manifests itself in the 
trilinear and quartic gauge boson self couplings obtained from Cq- 
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Leptons 



Quarks 



Families 
1 2 




T 




Ur 



R 



s'r 




R 



Qf 





-1 
-1 



+ 



+ 



T 



(3) 



+ 



+ 







w 



-1 
-1 
-2 



Table B.l: The classification of fermions into three families of weak isospin multiplets. 
The quantum numbers are the fermion charge Qf, the third component of the weak isospin 
7^ , and the weak hypercharge IV- The weak eigenstates {d',s',b') are obtained from 
the mass eigenstates [d, s, b) by Cabibbo-Kobayashi-Maskawa rotation. For each particle 
state, an antiparticle state which is obtained by charge conjugation exists. 



The Higgs Field and Higgs - Gauge Field Lagrangian 

It is a well-known fact that mass terms of the kind M'^B^B^ /2 violate gauge invariance 
(see e.g. ||8^) and are therefore forbidden. To have the carriers , W~ , and massive 
as is required by experimental observations, the SU{2)l x ?7(l)y,^ symmetry is sponta- 
neously broken via the Higgs mechanism [^], leaving behind an unbroken U{1) to be 
interpreted as the electromagnetic gauge subgroup U{l)em- The minimal Higgs mecha- 
nism introduces a complex scalar isodoublet field with = 1 



*(^) = % : (B.5) 




at the space-time point x. Its self-interaction and coupling to gauge bosons is described 
by 

£h = - (1?^$ ) - ) (B.6) 

with the covariant derivative 

V = d" - ig2TaW^'^'' + igl^B^' (B.7) 



and the Higgs potential 



V^($) = _^2|^$t<|,) +^($t$)' . (B.8) 
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The ground state of the quantum system with Higgs field is obtained by minimizing 
the potential V{^). In the ground state, the Higgs field has a non- vanishing vacuum 
expectation value which can be chosen as 



($)o = 





V 

V2 





2tM 



2A 



This reduces the manifest SU (2) symmetry to U{l)t 
as 



j • (B.9) 
The Higgs field can now be written 



V + H{x) + ix(x) 



/V2 



(B.IO) 



The fields 0^ and however, can be gauged away and therefore do not have any physical 
meaning. The remaining real part H{x) of 0° can thus be interpreted as a new scalar 
particle with its mass acquired via the quartic self coupling of the Higgs field: 

V2fi . 



M. 



H 



(B.ll) 



The kinetic part of the Higgs Lagrangian, eq. ( |B.6[ ), gives rise to mass terms for the gauge 
fields. Since Yw = 1 for the Higgs field, the gauge boson mass terms can be written as 







gi 











9l 



\ 




9i92 
9192 9l J 



(B.12) 



The physical parameters and fields are obtained from determination of the mass eigenval- 
ues and eigenvectors of equation ( |B.12|) . For the gauge boson mass part of Ch one thus 
finds 



^t7 



with the physical fields 

W^^ = 
= 

A^' = 

the masses 



v2 ^ 





M| 



cos^^ W^'^ + sinO^ 



(B.13) 



(B.14) 



w 



Mz 



2^2 



2V^i 




9l 



and the electroweak mixing angle defined by 

92 Mw 



cos9„ 



9i + 92 



(B.15) 



(B.16) 
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Identification of witli tlie plioton field coupling to the electric charge e yields an 
expression for e (compare equation (|B.14|) ): 

^^^^ fl'2sin6'^ = gicos6^ . (B.17) 



9i + 92 



The Fermion - Gauge Field Lagrangian 



As is seen from table fermions are grouped in left-handed doublets with field operators 




and right-handed singlets with field operators 

V'f = (B.19) 

with the family index j = 1, 2, 3 and the isospin index o" = ±. For the left-handed 
SU{2) doublets a is the sign of the third isospin component, and for the right-handed 
U(l) singlets it denotes the corresponding particle type. The color index of quarks is 
suppressed. The Lagrangian for the fermion - gauge field interaction can be written in 
terms of the covariant derivative (see equation ( [B.7| )) 



with 



The latter definition ensures that if) is identical to the commonly used one. 
The Yukawa Interaction 

Because of the different gauge transformation behavior of left- and right-handed fermion 
fields, mass terms of the kind mipil) are forbidden. Fermion masses can, however, be 
dynamically generated via the gauge invariant Yukawa Interaction 

Cvuka^a = - E {[^?i-#'^^f- + g.+ i^f+^^n^"^] + h.C. \ . (B.21) 

j 

Here, we have used the Pauli matrix 



The fermion masses are derived from the Yukawa coupling constants gja and the non- 
vanishing vacuum expectation value of the Higgs field operator. The Yukawa Interaction 
takes on a very simple form in the unitary gauge where the unphysical degrees of freedom 



of the Higgs field are gauged away 



^Yukawa = 'Y.^f'^f^f " E ^f'^f^ ' (B.23) 
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where the sum runs over all fermions. Note, however, that this is a description in a specific 
gauge and not a gauge-invariant Lagrangian density. 



Quantization of the Lagrangian, Gauge Fixing and Ghosts 

In general, to be able to quantize the Lagrangian described above, it is necessary to apply 
a gauge fixing procedure. Simultaneously, a compensation of the unphysical content 
of the gauge fixing Lagrangian Cqf is needed. This compensation is provided by the 
introduction of a Faddeev-Popov ghost term Cpp- Details of the procedure can be found 
in references 



26, 86 



or in standard textbooks |8^. As this thesis is worked out in 
the unitary gauge, no additional gauge fixing and ghost terms are needed, no unphysical 
particles appear in the actual calculation. As only QED corrections are treated, the 
subtleties connected with unitarity and renormalizability of the unitary gauge do not 
play any role. 



B.2 Feynman Rules 

One can express the Lagrangian as exposed in appendix [B.l| in terms of physical para- 



meters and derive the Feynman rules for the Standard Model. Below, Feynman rules 
relevant for this thesis are presented. Complete lists of Feynman rules are found in the 
literature [^, |8^, Q and in the unconventional, but excellent textbook by Veltman ||91[ . 



Rule 1: External Lines 

External fermions with momentum p and helicity index r, external gauge bosons with 
momentum k and polarization index A. 

incoming fermion : u'^{p,m) 

outgoing fermion : u'^{p,m) 

incoming antifermion : 'ir{—p,m) 

outgoing ant if ermion : u''[—p,m) 



p, 


r 


p, 


r 


p, 


r 


p, 


r 



a 

a 



mcommg 



outgoing 



gauge boson: 6 '"(k) 



The fermion spin and vector boson polarization summation conventions are given by 

r 

r k^ 



J2£^''"{k)e^'^{k) 

A 



(B.24) 
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Rule 2: Propagators 

Propagators with momenta p for fermions, k for photons, and q for bosons. The imaginary 
parts —ie ensure proper causal behavior. Define = ^y^Pfj,. 



fermion 



+ im 



p p^ + — i£ 

Q; P i q^g 

•J\jxr\j^u\j\j\j\j\j* photon : — — 

k k^ — IS 

« • K V i + ?a W^V) 

• * • massive gauge boson V : - 



+ - \£ 



Rule 3: Vertices 



In this thesis, only neutral current vertices are needed. They can be cast into a simple 
generic form: 



r 



B ^-f, -93 92- T {g]fB + afjB 75) 



For the couplings 9b-i dJfB^ dfjB fermion / to the photon and the the 

relations 

e 

sin^^ 

9i ^ ^ 9z ^ ^ 



4 ■ cos6l. 



w 



9ff^ = sin6'^ Qf 9^^^ = 2T^^} - 4 • sin26'^ Qf 

hold, with the fermion's electric charge Qf in units of the positron charge e and the 

(s) 

fermion's weak isospin third component T^r f ■ It is often desirable to write the vertices 
in terms of left- and right-handed couplings, 

- 9B92- r {gm + 9mi^) = - 7" + 75) + - 75)) (B.26) 

with 

LffB = (g]jB + gfjB) RffB = (g^fj^ - 9fj^ . (B.27) 
For the above neutral current vertices one obtains 

= 2 = 2 sine! cosO, (4!}-sm'»w*) 

^'A = ^ "''^^ - 2 sino! cose. ^'"'"•"^ ■ <^-^«) 
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Rule 4: Ordering 

The spinor factors (7-matrices, propagators, 4-spinors) for each fermion hne are ordered 
from right to left while following the line in its arrow sense. Multiply the expression by 
a phase factor +1 (-1), if an even (odd) permutation is required to write the fermion 
operators in the correct normal order. 



Rule 5: Four-Momentum Conservation 

The four-momenta of lines meeting at a vertex satisfy four-momentum conservation. For 
each four-momentum q which is not fixed by four-momentum conservation, carry out the 
integration J d'^q. 



Rule 6: Fermion Loops 

For each closed fermion loop take the trace and multiply by -1. 



Rule 7: Symmetrization 

Introduce a relative sign +1 (-1) for diagrams that are obtained by an even (odd) number 
of permutations of identical final state fermions. Obtain the matrix element by adding 
all diagrams. 



Rule 8: Differential 2 — »n Cross-Section 

Aiif be the matrix element for the scattering of a two-particle initial state i = {1,2} to 
an n-particle final state / = {3, 4, ...,n + 2}. Then obtain the differential cross-section by 



\M 



da = {2ny 5^'\P, - Pj) dp, ■ ■ ■ dp^^, 

A J (pi-p2y — 77111712 ^ 



n+2 

Pi = P1+P2 Pf = Y.Pk ■ (B.29) 

fe=3 

S yields the symmetry factor, if there are I sets with ki identical particles, 

s = U^ii . 

I 

With p^ representing the particle energy, p^ = y/p'^ + m^, the Lorentz invariant phase 
space element is given by 

d^p 

dp - 



(27r)3 2/ 



In case of unobserved polarization degrees of freedom, the average over initial state and 
the sum over final state polarizations has to be taken. The factors 2m left out in equa- 
tion ( [B.29| ) compared to common textbooks like |^ is cancelled by the normalization 
chosen for the fermion spinors u{p,7n). 
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B.3 Renormalization 



In higher order perturbation theory, the relations between the parameters of the Stan- 
dard Model Lagrangian ( [B.2D and the physical observables differ from tree level. At loop 
level, the "bare" parameters of the tree level Lagrangian have no more physical mean- 
ing, and a redefinition of parameters and fields is required to recover physical meaning. 
The procedure of redefinition is called renormalization and properly removes divergences 
appearing in loop diagrams. For a mathematically consistent treatment of divergences, 
the theory must be regularized which is most commonly done with dimensional regular- 
ization. The "renormalizability" of the Standard Model assures that a finite number of 
so-called renormalization constants suffices to redefine physical quantities in all orders of 
perturbation theory Renormalization schemes differ in the choice of input parame- 



ters and renormalization conditions. As physical results must be scheme independent, the 
observables calculated in different renormalization schemes are identical in infinite order 
perturbation theory. Results from ra*'* order perturbation theory obtained in different 
schemes, deviate from each other in higher order contributions |^|. The accuracy of the 



finite order approximation is influenced by the choice of input parameters. Among the 
various renormalization schemes found in the literature, the most prominent ones are 

• The "modified minimal subtraction" (MS) scheme [p3| , where the role of the 



electroweak mixing angle is emphasized. Renormalization constants are fixed by a 
simple subtraction of the singular parts 

^ + ^ - ln(2 V^) 



n 



from the two- and three-point functions calculated in dimensional regularization (see 
appendix ^ . The common input parameters for the MS scheme are 

a, Mw, Mz, Mh, and the fermion masses mj. 

An advantage of the MS scheme is that theoretical uncertainties due to quark in- 
duced self-energy effects are much reduced p4 |. 



The * scheme which emphasizes vector boson propagator effects and running 



parameters |^6[. The * scheme is not very commonly used today 



The Gu scheme |92, M uses the input parameters 



G^, Mw, Mz, Mfj, and the fermion masses nif. 

It is important that does not run from low energies up to the vector boson mass 
scale. The fine structure constant a and the Fermi couphng constant G^ are related 
via Ar, the non-QED correction to /i-decay. 



'''In "unrcnormalizable" theories, additional renormalization constants are needed for each new per- 
turbation order. 
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• The on-shell scheme p6| , |98| which was used for this thesis and will now be described 
in some detail. 

In the on-shell renormalization scheme, to recover physical meaning at loop level, so-called 
multiplicative renormalization of bare fields 0o aiid bare couplings Qq is applied, 

00 = ^70 = Zgg . (B.30) 

The renormalized parameters g are finite and fixed by renormalization conditions. Field 
renormalization ensures finite Green functions |^. Using the expansions Zi = l + 6Zi, 
the bare Lagrangian is split into a renormalized part and a counterterm part, 

£(00, go) = C{Z'J'4>, Zgg) = £(0, g) + 6C{(P, g, 6Z^, 6g) . (B.31) 

The renormalized Lagrangian can be expressed in terms of the chosen set of physical 
parameters, and from 6C one derives counterterm Feynman rules in accordance with the 
on-shell renormalization conditions. Renormalization conditions fix the renormalization 
constants and thus give meaning to the theory. Counterterms properly remove ultravio- 
let divergences arising from loop diagrams. Infrared divergences from loop diagrams are 
cancelled by real bremsstrahlung radiative corrections. In the on-shell scheme, renormali- 
zation conditions are chosen to have propagator poles at the physical particle masses with 
residue = 1 and to recover the classical Thomson limit for the ee7-vertex. The common 
set of input parameters for the on-shell renormalization scheme is 

e, Mw, Mz, Mh, and the fermion masses m/. (B.32) 

It is a vital property of the on-shell scheme that tree level relations between masses, 
couplings and the electroweak mixing angle 6w remain intact also at loop level. Thus one 
may easily switch between input parameters. In particular 

sin^^^ = 1-^ (B.33) 

can be retained in the on-shell renormalization. For the renormalized Fermi coupling 
constant which is an attractive input parameter, because it is experimentally known 



with high precision, the following relation holds |96 



Thus, in the case of pure QED loop corrections to a process, even retains its tree level 
relation to the other parameters. 

As this thesis is concerned with QED corrections, the on-shell renormalization for QED 
will be worked out in some detail below. After choosing a set of input parameters, one 
starts from the tree level QED Lagrangian 

F^"" = d^'A^Q-d'^A^ (B.35) 
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with the photon field Aq and introduces multiphcative renormahzation constants: 
V^o = Jz^-^ = Jz'a-A^ 



mo 



■ m 



eo 



Therefore, expanding Zi = l-\-6Zi, the QED Lagrangian is given by 
with the counterterm Lagrangian 



6C 



(B.36) 



(B.37) 



6C 



4 



(B.38) 



From the requirement of gauge invariance one concludes that the counterterms SZ^ and 
6ZjIj are equal which is nothing but the QED Ward identity in the language of renor- 
mahzation constants |88]. In the next step, counterterm Feynman rules are derived from 
6C. Counterterms are then added to physical 1-loop amplitudes and finally fixed by the 
on-shell renormahzation conditions. The QED counterterm Feynman rules are given on 



page 50 after explicit expressions for the counterterms have been derived. 



The photon field counterterm SZa 

Renormahzation of the photon propagator means to add an ultraviolet divergent vacuum 
polarization diagram and a counterterm diagram as given in figure |B.1| to the photon 
propagator derived from the renormalized Lagrangian C in equation ( p.37|) . From the 
three diagrams in figure |B.1| , the 1-loop expression for the photon propagator is derived 
as 



kP' — is 



+ 



-1^ 



— ie 



in 



a/3 



-1^ 



I3u 



k'^ — ie 



(B.39) 



with the vacuum polarization given by (see e.g. p9| ) 



(27r)" 



Tr 



7o 



+ im) 7/3 - ^ + i"^) 



+ m? — ie 



ip — ky + — ie 



fj.=me. 



4i (eQf)^ 
16 7r2 

iga^Wie) + ik^kpU'^'\e) 



2 

o [k'^gap-kakp) P + 2 {k^gai3-kakfs) I{k'^,m^,m'^) 



(B.40) 



-7-/ 7 2 2 2 

l{k , m , m 



dx x(l — -In 



m 



x(l -x)— H 



(B.41) 
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a) 



p,m 




(p-k),m 
b) 



c) 



Figure B.l: Feynman diagrams for the photon propagator renormalization. a) Renormal- 
ized photon propagator, b ) Vacuum polarization, and c ) Photon propagator counterterm. 



The dimensionally regularized ultraviolet divergence P is defined in equation ( p.7| ). It 
follows from current conservation that, in equation ([B.40I ), the term proportional to the 
photon momentum k yields vanishing contributions p8l. Therefore one finds 



-ig 



A;2 



1 + 



W - kHZA 



(B.42) 



Introducing the renormalized vacuum polarization fl''' 
quantity, one can rewrite equation (p.42| ): 



W - k'^ 5Za which is a "small" 



k"^ — ie 1 _ 

-ig'"' 
k'^-ie- t\r< 



+ higher orders 



higher orders 



(B.43) 



It is seen from equation (|B.43|) that the photon propagator has a pole strictly at fc^ = 0, 
if TV = k'^d'W / dk'^\k'i=Q +0{k^). Then, the photon does not acquire a mass term through 
higher order effects and hence there is no photon mass renormalization. The on-shell 
renormalization condition requires the photon propagator (|B.42|) to have residue = 1, i.e 
to retain its tree level form in the limit k"^ = 0. Thus U'' /k"^ 



fc2=0 



: and therefore 



5Za 



k^ 



(B.44) 



fc2=0 



which is readily evaluated from equation (|B.41|) , yielding 



[eQff 



.1 2P 



mi , 



(B.45) 



167r2 3 

For non- vanishing fc^, fl''' yields the charge renormalization familiar from many text- 



books 107, 



rog. 
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(k-p),m 

k,m 1, k,m k,m 

• » • + • » iAy\y\yv» » • + • )( • 

p 

a) b) c) 

Figure B.2: Feynman diagrams for the fermion propagator renormalization. a) Renormal- 
ized fermion propagator, b ) Fermion self energy, and c ) Fermion propagator counterterm. 



The fermion field and mass counterterms 6Z^, SZm, and 6m. 

The Feynman diagrams for the renormahzed fermion propagator are shown in figure |B.2 . 
Decomposing the fermionic terms from the counterterm Lagrangian, 



im 6Zr, 



im) 6Z^ — i Sm 



5m = m {6Zm — 5Z^) , (B.46) 
the 1-loop expression for the fermion propagator corresponding to the three diagrams is 



-im — le 



i 6m 



1 



im 



-im — ie 



^ —im—ie 



(B.47) 



with the fermion self energy given by 



- i ^ im) 7^ 



d^p 

WT T(A;-p)2 + m2 - 



S^(im) + (,2;; - im) + O [(^ - 



im 



S'^(im) 



m 



16 7r2 

{eQfY 
16 vr^ 



im 6P - 4 + 3 In — 



mi 



2P + 4P 



IR 



4 + 3 In ■ 



m 



mf 



(B.48) 



The dimensionally regularized infrared singularity P^^ is defined in equation ( |G.8| ) . Using 
the fermion self energy decomposition from equation (|B.48|) , the renormalized fermion 
propagator can be written as 



— im — 



1-6Z^ + B^ + 0{I^ - im) + 
1 



i 6m + S-^(im) 



— im 



(B.49) 



(^-im) 1 + 6Z^-Bf - 0{l^ -im) -i(5m-S/(im) 



+ higher orders . 

(B.50) 
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Imposing the on-shell condition, i.e. requiring in equation ( |B.50D to have a pole at 
^ = im, the result for the fermion mass counterterm reads 



6m 



im 



m 



16 7r2 



6P 



, m 
4 + 3 In — 

mi , 



(B.51) 



In addition, as the pole of the 1-loop fermion propagator is required to have residue = 1, 
the second fermion propagator renormalization condition reads 



16 7r2 



2P + 4P 



IR 




(B.52) 



With the renormalization conditions ( |B.51| ) and ( [B.52| ) it is ensured that external on-shell 
fermion lines in Feynman diagrams do not receive any QED loop corrections. The mass 
counterterm dZ^ is now easily derived by introducing the results ( p.51|) and (p.52| ) into 
equation ( |B.46| ): 



5Z„ 



6Z,i 



6m 
m 



16 7r2 



8P + 4P^^ - 8 + 6 In 



m 



mt 



(B.53) 



The vertex counterterm 6Ze 

The renormalized 1-loop electromagnetic vertex consists of the three Feynman diagrams 
depicted in figure p.3| . For on-shell external fermions it is given by 

- eQ/f^ = -eg/7^(l+5Ze) - eQ/F^ (B.54) 

with the vertex correction 

-n) [ ^> 7''(^2-|*+i"^)7^(^i-^+i"^)7M 



F^ = -i{eQff f.^'--^ J- 



16 7r2 ^ ' 



(27r)" (^/^^_p)2_|_^2_j^ p)2-|-m2— ie 

2P - 2 (q^ + 2m') Jiq\ m\ m^) P^^ - In ^ + 



p^ — ie 



+ ^(7^ -^7^) J{q',m',m^) 



^^■Fi{q^,m^,m^) + (7V " ^7^) -^2(9 



2 2 2n 

,m ,m ) 



(B.55) 



where 



dx 



2=n 1 



q^=0 



x{l — x)q'^ + 777.2 



777^ 



dx 



x{l — x)g2 + rn? 



xil — x)q^ + m^ q2=Q 1 7772 
■ In = — - ■ In 



777 f 



7772 7772 



(B.56) 
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a) b) c) 

Figure B. 3: Feynman diagrams for the electromagnetic vertex renormalization. a) Vertex, 
b) Vertex correction, and c) Vertex counterterm. 



From equation ( p.55|) , the Dirac form factor Fi and the Pauh form factor F2 can easily be 
read off. The on-shell renormahzation condition requires the renormahzed vertex F^ to 
recover 7^ in the Thomson hmit g — > 0, 

f^(g = 0,m^m2) = 7^ . (B.57) 

Equivalently one may require 

= o,m^m2) = -5Ze (B.58) 

for the Dirac form factor. Introducing condition ( p.57|) or ( [B.58|) into equation (p.54| ) and 



using the result for F from equation ( pj.55| ) with the proper limits J{q = 0, m , m ) and 
/C(g^ = 0, m?, m?), one obtains 

5Ze-7' = -F\0, m2,m2) = Ml! f 2P + 4pi^ - 4 + 3 In ^ V 7' (B.59) 

for the vertex counterterm. It is therefore seen that 5Z^ and 5Zf. are equal and thus 
indeed fulfill the QED Ward identity. 

It should be mentioned that the vertex counterterm derived above is also valid for the QED 
renormalized ffZ^ vertex. Briefly, one could argue that the term (yO^fz + dffz ^s) fro^ 
the replacement 

- eQ/ 7^ — ^ - 9z92- 7^ [a^fz + ajfz Ts) 
is exactly pulled through to the end of the vertex correction, 

= -9z92-T^-{9jfz + 9ffzl5) = 4 ging^'cosg^ ' ' {dffz + dffz I5) • 

(B.60) 

Therefore, as the weak mixing angle is not affected by the on-shell QED renormalization. 



49 



Thus the renormahzation condition ( |B.57| ) and its solution ( |B.59| ) yield the same multi- 
plicative counterterm for the ffZ^ and the ffj vertex. 



The same result may be obtained from equation (58) of reference |[101|| or from equation 
(4.21) of reference ||102|| by neglecting the counterterms that do not belong to the on-shell 
QED renormahzation. 



Counterterm Feynman rules 

With the above results, namely equations ( |B.45D , ( |B.51[ ), ( [B.52| ), and ( |B.59D , the coun- 
terterms as needed in the 1-loop QED counterterm Feynman rules are fixed. The QED 
counterterm Feynman rules are listed below and must be used together with the set of 
rules given in section [B.2| . Subsequently the abbreviation c.t. will be used for "counter- 
term" . 



Rule 9: Renormalized external lines 

Do not add loop contributions nor counterterms to on-shell external particle lines. 



Rule 10: QED Counterterms 



P 



fermion c.t. : 



-im — ie 



■im) 6Z^ — iSm 



im — ie 



a P 
k 



/■ 



a 



photon c.t. : 



P-i 



k'^ — ie 



k'^ — ie 
- 5Za\ 



7, Z° vertex c.t. : -9392-1°' {9jfB + 9ffB Ts) ■ SZe 



From the above counterterm Feynman rules, several observations are made: 

• In case of electron lines, the logarithms \n{m?/m?^) vanish in all counterterms. 
Initially, these logarithms arose due to the choice of the electron mass as the scale 
for dimensional regularization. 

• For massless fermions, the fermion propagator counterterm Feynman rule translates 
into a factor —5Z^ that must be multiplied to the tree level amplitude corresponding 
to the counterterm amplitude under consideration. 
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• Similarly, an amplitude with a photon propagator counterterm is obtained by mul- 
tiplication of the corresponding tree level amplitude with —SZj^. 

• An amplitude with a QED vertex counterterm is obtained by multiplication of the 
corresponding tree level amplitude with SZ^. 

In conclusion, the counterterms for the QED initial state corrections to process ( |1 . 1| ) may 
be multiplicatively taken into account as it is done in equations ( pD.lSj ) and ([E.16|) . 
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Appendix C 

Complex Logarithm and 
Polylogarithms 



For the proper computation of many phase space and loop integrals, the definitions of 
the complex logarithm and the derived polylogarithms are needed. In this appendix, 
the necessary definitions are listed together with some useful relations. An exhaustive 



treatment of Polylogarithms is found in reference ||103 



C.l The Logarithm Function 

Throughout this dissertation the logarithm is defined as the main branch of the inverse 
of the complex exponential function: 



Inz 



ln(| 



z\ ■ e 



iarg(2) 



= Inl^l + iarg(^) , z e C\iR~ 

aig{z) e ]-TT,n[ (C.l) 

The logarithm of a product of two complex numbers a and b is given by ||104| 

ln(a ■ 6) = In a + In 6 + ri{a,b) 



(C.2) 



r]{a,b) 



2m 



G(— a) G(— 5v»Ti b) Q(^sm ab) — B(Qm a) Q{^ism b) 0(— ^vsn ab) 
with the Heavyside step function Q{x). It can be concluded that 

ln(a ■ 6) = In a + Info if $)ma and '^mb have different sign. 

In- = In a — In 6 if ^>ri a and '^rn b have equal sign. 





From table |C.1| , where all possible sign combinations for ri{a, b) are listed, it can be seen 
that for real numbers x and y and infinitesimal e 



ln(x?/ — ie) 



ln(x — ie) + ln(?/ — ie/x) . 



(C.3) 
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sgn(5v5ma) 


sgn(9m6) 


sgn(^>n ab) 


r]{a,b) 


+ 


+ 


+ 





+ 


+ 




-2ni 


+ 




± 







+ 


± 









+ 


+27ri 












Table C.l: All possible values ofri{a,b). 



C.2 The Dilogarithm Function 

The dilogarithm function is defined by the integral 



ln(l - tz) 
t 



dt 



z e€\{xeR: X > 1} . 



(C.4) 



The following one parameter relations are valid for z being any complex number such 
that none of the involved arguments lie on the corresponding logarithmic or dilogarithmic 
cuts. 

Li2(-ll 



Li2(0 

U2{Z 
U2{Z 
U2{Z 
U2(Z 



= 0, 

+ U2{-Z) 
+ Li2(l/^) 

+ Li2(l-2; 
+ Li2 



TC 

12 



Li2(l) = y 



TT 

y 

7^ 

~6 



In' 



\n{z) ■ ln(l - z) 



1 



l-l/zj 
Li2(l-z) + Li2(l-l/2) 



- In^ri-^) 
2 ^ ' 



1 



In'(z) 



(C.5) 
(C.6) 
(C.7) 
(C.8) 
(C.9) 
(C.IO) 

(C.ll) 



For an infinitesimal real number e and xeM the dilogarithm satisfies 

{Li2(x + ie)} = sgn(e) ■ 9(a; — 1) • vr ln|a;| . 

Therefore, for xeM, x > 1 one can easily deduce 

Li2(a; + ie) + Li2(x — ie) = sgn(£:) ■ 27ri ■ Inx . 

A general two parameter relation for complex- valued dilogarithms which can be useful for 
special purposes and will not be repeated here is given in reference ^6 . 
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C.3 The Trilogarithm Function 

Similar to the dilogarithm, the trilogarithm function is defined by an integral: 



dz 



hiz ■ ln(l — yz) 



dz 



y e €\{xeR: x > 1} . 
From this definition two integrals are easily derived: 



(C.12) 



In z ■ ln(l — z) 



dz = U-i{y) - \ny-U2{y) 



ye€\R U [0,1] (C.13) 



In 2; ■ ln(l + az) 



dz = Usi^-ay) - \ny ■U2(-ay) , 

ye€\R-, ay e €\{x e R : x < -1} 



(C.14) 



Some useful one parameter relations are given below. Again, it is understood that the 
arguments in each of all equations must not lie on any of the corresponding above-defined 
cuts: 



Li3(0) = 

Li3(l) = C(3) = 1.20205690... 

1 



Li3(-1) = --Li3(i: 



U3(z) + LU-z) 



U,iz') 



(C.15) 
(C.16) 

(C.17) 



Us{z)-Us{l/z) 



IT 



ln(-;z) 



6 



In^f-^) 



(C.18) 



Further, more complicated relations are found in chapter VI of reference ||103|| . Here, only 
some relations for real arguments x will be quoted in addition: 



TT^ 1 1 

LU(x) — LU(l/x) = — Inx — - In^x — -ivrln^x 



X > 1 (C.19) 



Li3 



-X 



1 — X 



+ Li3(l -x)+ Li3(x) 



Li3(l) + — In(l-x) - - Inx-ln^(l-x) + -ln^(l-x) , 0<x<l (C.20) 
6 2 6 



Lis 



-X 



1 — X 



Lis 



1 — X 



Li, fx) 



Li3(l) - y ln(l -x) - ^ In(-x) ■ln^(l - x) + ^ ln^(l - x) , x<0 . (C.21) 
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Appendix D 

Description of the Phase Space 



It is the intention of this appendix to famiharize the reader with the phase space pa- 
rametrizations used throughout this thesis. For the Born case of process and for 

the virtual initial state corrections, the 2 — > 4 particle phase space is needed, whereas the 
2 — » 5 particle phase space parametrization is used to describe bremsstrahlung. 



Figure D.l: Graphical representation of a 2^4 particle reaction. In general, the vectors 
Pi,P2,P3 and Pi do not lie in the plane of the picture. 
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D.l The 2^4 Particle Phase Space 

A reaction of two particles to a four-particle final state has eight kinematical degrees of 
freedom which can be parametrized in many different ways. For semi-analytical calcula- 
tions the following choice of kinematical variables has proven to be convenient |6^ : 

: Azimuth angle around the e"*" beam direction k2- 

: Scattering angle of Vi with respect to the e"*" direction k2 in the center of 
mass frame. 

Azimuth angle around the Vi direction. 

Decay polar angle of pi in the Vi rest frame with the z axis along Vi. 
734 . Azimuth angle around the V2 direction. 
^34 : Decay polar angle of in the V2 rest frame with the z axis along V2. 
S12 : Invariant mass of the final state fermion pair /i/i : S12 = —vf = —{pi +P2)'^- 
S34 : Invariant mass of the final state fermion pair 72/2 : S34 = — f| = —{ps +^4)^. 

Here, ki and k2 are the momentum four-vectors of the initial state electron and positron, 
Vi and V2 denote the four-momenta of the intermediate bosons and Pi,P2,P3, and p4 are 
the final state four-momenta. The corresponding three-momenta are given by superscript 
arrows. To illustrate the above kinematical variables, a graphical representation of a 
2 — »• 4 particle reaction is shown in figure p.l| . The center of mass energy squared is given 



(Pl2 
O12 



by s = — (fci + A;2)^ = — (pi + P2 + Ps + ^4)^- The decomposition of the 2 — > phase 
space into subsequently decaying compounds is, for example, worked out in chapter 4.2 of 
reference [ 105|| . In terms of the above variables the four-particle Lorentz invariant phase 



space has the form 

,-p 1 d^Pi d^P2 d^P3 d^pi ^ 



1 y^A(s, 512,534) ^A(si2,mi,mi) \J \{s-n,ml,ml) 



(27r)ii 8s 8si2 8S34 

dsi2 ds^i dcosd d(f)i2dcos9i2 d(j)^4,dcos9^i , (D-1) 
with the usual definition of the A function, 

A(a, 6, c) = + l)^ + - 2ah - 2ac - 2hc 
which is symmetric in all arguments and has the special cases 
\{a,h,h) = -(1-46/0) , 
A(a,6,0) = {a-bf 
A(a,0,0) = 
The limits of the phase space variables are: 





< 


S12 


< 




-2m2)' 


Ami 


< 


■S34 


< 




- V^u) 


-1 


< 


cos'd 


< 




+1 





< 


0{12,34} 


< 




27r 


-1 


< 


C0S6'{12,34} 


< 




+1 
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In the presented semi-analytical cross-section calculation, the integrations over the final 
state fermion decay angles 0i2, 6*12 and 034, 6^34 are separated from the other integra- 
tions and are carried out analytically with the help of invariant tensor integration (see 
appendix |G.1| ). The remaining integrations over {}, S12, and S34 are non-trivial, especially, 
if virtual QED corrections are included. The angle is integrated analytically, S12 and 
S34 are integrated numerically. To perform integrations over the above variables, one 
must express the four-vectors of the four-particle final state in terms of the phase space 
parameters. All scalar products appearing in the cross-section calculation can then be 
expressed in terms of these integration variables. In the center of mass frame with a 
Cartesian coordinate system as drawn in figure 
given by 



the initial state vectors ki and ^2 are 



ki = {k^, k sin'd, 0, — k cos-i?) 
k2 = — k sint?, 0, k cost?) 



with 



k 



2 



(D.2) 



4m2 



and the electron mass rrie- Using the mass mi for /i, and m2 for /2 one obtains further 



Pi = {luPif + I12P12 cos6'i2, P12 sin6'i2 cos0i2, P12 sin6'i2 sin0i2, 7i2Pf2 cos6'i2 + 7i2pf2° 



P-i = (712^2° - 7i2pf2 cos6'i2, -Pi2 sin6'i2 cos(/)i2, -Pi2 sin6'i2 sin0i2. 



B,0 



I12P12 cos6'i2 + 712P12 



P3 = (734^34'° - 734pf4 C0S6'34, pf^ Smd^^ COS(/)34, ^^4 sin6'34 Sin(^34 , 734^9^4 cos6'34 - 



B,0 

734P34 



Pi = {luPsf + IsaPu C0S^34, 



-p^4 sin6'34 COS034, -pgt sin6'34 sin034, 

-luPii COS^34 



B,0 
iMPzi 



with 



(D.3) 



Bfl 

Pli 



Bfl 

P34 



P?2 



7?2 



7l2 



^J\{si2,ml,ml) 

S + Si2- S34 



'2y/SSi2 
\J S12, S34) 



pI 



73°4 



734 



^A(s34,mi,mi) 



12 



2a/S34 
S — Si2 + S34 

A(s, S12, S34) 

2a/SS34 



The superscript B hints at the relation of the corresponding quantities to the appropriate 



boson rest frame. The quantities 



P 



B 



, and 7°- are e.g. found in chapter IV, 



equations (3.4) and (3.6) of reference ||106|| . 7ij = |7ij| may be derived from the equation 
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7°/— 7i7- = 1- use in appendix p.2| vector components in the above equation ( p.3| ) are 



It] 

rewritten in terms of the generic names 



Pt = P-is), pKs), ptis)) . (D.4) 

All possible scalar products built from the four- vectors ki, k2,Pi,P2,P3, and p^ can now be 
easily expressed in terms of the phase space parameters and the particle masses. Scalar 
products involving Vi or f 2 are derived from Vi = pi + p2 and f 2 = Ps + Pa,- 



I 



V2 



S + S12 — S34 



, 0, 0, 



A(s, S12, S34 



V 



2yi 



2yi 



(D.5) 



Subsequently, the fermion masses mi and m2 are neglected. The initial state electron mass 
rUe is neglected wherever reasonable which is called the "ultrarelativistic approximation" 
(URA). The following quantities are often used throughout the calculation: 

A = A(s, 312,534) , (T = S12 + S34 , 6 = S12 — S34 , 



- s — cr 



= i (s - a + 



(D.6) 



D.2 The 2^5 Particle Phase Space 

In the case of initial state photon bremsstrahlung a fifth particle appears in the final state. 
The five-particle phase space has eleven kinematical degrees of freedom which are chosen 
as 

: Azimuth angle around the e"*" beam direction ^2- 

— * 

9 : Polar (scattering) angle of the photon with respect to the e"*" direction k2 in 
the center of mass frame. 
(f)R : Azimuth angle around the photon direction p. 

Or : Polar angle of Vi in the {V1+V2) rest frame with the z axis along {vi+V2)- 
012 : Azimuth angle around the Vi direction. 

^^12 : Decay polar angle of pi in the Vi rest frame with the z axis along Vi. 
034 : Azimuth angle around the V2 direction. 

^34 : Decay polar angle of p^ in the V2 rest frame with the z axis along V2. 
S12 '■ Invariant mass of the final state fermion pair /i/i : S12 = —vf = —{pi +p2)'^- 
S34 : Invariant mass of the final state fermion pair 72/2 : S34 = —V2 = — (ps +^4)^. 
s' : Reduced center of mass energy squared. This is equivalent to the invariant 
mass of the final state four-fermion system: 
s' = -{vi + V2y = -{P1+P2+P3+ Pif- 
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Figure D.2: Graphical representation of a 2-^5 particle reaction. In general, the vectors 
Vi,V2,'Pi,'P2,'P3, and Pi do not lie in the plane of the picture. In the figure, the angle 9 is 
called 9^ to emphasize that it belongs to the direction of the radiated photon. 



The symbol p denotes the four-momentum of the photon radiated from the initial state. 
The kinematics of a 2 — 5 particle reaction is graphically illustrated in figure p.2| . The 
center of mass energy squared is given by s = — (fci + /c2)^ = —{p + Pi+Pi+Pz+ViY'- The 
angle 9r corresponds to the angle t9 in the 2^4 phase space described in appendix |D.1| , 
where the two boson rest frame and the center of mass frame of the e"'"e~ annihilation are 
identical. Due to the radiated photon, the two boson rest frame and the e"'"e~ center of 
mass frame do not coincide for the 2—^5 phase space. The five-particle Lorentz invariant 
phase space is parametrized by a decomposition into subsequently decaying compounds, 

_ 1 s- s' y^A(s',Si2,S34) y^A(si2, mf,mf) ^\{s:n,ml,ml) ^ 
^ ~ (27r)i4 ~8s 8s' 8s^2 8s^4 ^ 

ds' dsi2 dss4, dcos9 d(j)ji dcos9R d(j)i2 dco?,9i2 dcf)^^ dcos9u (D.7) 

with phase space limits 
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(2mi + 2m2) 



Ami 

Ami 
-1 

-1 





r,' 

s 




S 


< 


Sl2 


■ 


( r~i r> 

\^Js' - 2 


< 




< 


(v^-v 


< 


cos6 


< 


+1 


< 


0{R,12,34} 


< 


27r 


< 


C0S6'|^_12,34} 


< 


+1 



12 



(D.8) 



or alternatively 



Am\ 
Am\ 



-1 

-1 



< 


S\2 


< 




- 2 777,2)^ 


< 




< 






< 


s' 


< 




s 


< 


cos6 


< 




+1 


< 


0{iJ,12,34} 


< 




27r 


< 


C0S6'{/J^12,34} 


< 




+ 1 



(D.9) 



In the semi-analytical bremsstrahlung cross-section calculation, after integrating over 
012, 612 and 034, ^34 with invariant tensor integration, the non-trivial angular integrations 
over (pR, 6ji, and 6 are performed analytically. The invariant masses S12, S34, and s' are 
then integrated numerically, using the parametrization ( p.9| ). Below, the four- vectors of 
the five-particle final state are expressed in terms of the variables of the phase space ( p.7|) . 
Vectors are given for the center of mass frame with the Cartesian coordinate system of 
figure p.2| . Similar to equation (p.2| ), ki and ^2 are represented by 



ki = (A;°, k sin9, 0, — k cos9) 
k2 = {k^, — k sm9, 0, k cos9) 



with 



k 



2 
2 



(D.IO) 



= a/1- 



4m,^ 



Secondly, one finds 
0, 0, 



P 



with 



P 



2v^ 



(D.ll) 



for the four-momentum of the photon radiated along the z axis. The momenta of the 
final state fermions can be written down in terms of the expressions originally introduced 
for the four-particle final state in equation (p.4|): 



Pi 



1 



(0) 



\ 



pI{s') sm9R + pf{s') cos6'Rj coscpR + p\{s') smcpR 
pl{s') sin6'H + Pi{s') cos6'r) sin^^j + p\{s') cos(j)R 
y 7(0) (p^^{s') cos9r - (s') sm9R^ - 7P°(s') 



(D.12) 



with 



7 



(0) 



s + s' 



2Vss' 



7 



2Vss' 



(D.13) 
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In equation ( p.l2| ) the components of the momentum four-vector have been written in a 
column for reasons of clarity. Equation ( p.l2|) is better understood, if one notes that, in 
the Pi + P2 + Ps + = Vi + V2 rest frame, the fermion momenta are given by 



Pi 



P-is'), pKs'), Ptis'] 



(D.14) 



where, compared to equation (p.4| ), only the center of mass energy had to be adjusted 
due to photon radiation. After rotating and boosting to the five-particle center of mass 
frame, equation ( p.l2| ) is obtained. 

The intermediate boson momenta vi and V2 are given by expressions similar to those in 
equation (p.3|) , 



V2 



^(0) pR^o _ cos9r, Vr smdR cos0r, - Vr smdR sin^R, 

'J^^^VrCOsOr - 'jpf.f 

(7^°^ P34° + I'Pr cosOr, - Vr sinOR cos0ij, Vr sinOR sin0R, 

-7(°)pRCOs^fi-7P34°) . (D.15) 



with 



V, 



•^(s'; ■Sl2; -534) 
S' + S12 - S34 



2Vs 



R,0 



S - S12 + S34 



(D.16) 



The four-momentum of the two-boson system recoiling against the radiated photon is 
given by 



{V1 + V2) = (7^°W, 0, 0, -7V^) 



0, 0, ^1 



(D.17) 



Once more, some derived kinematical quantities are defined for use throughout the cal- 
culation: 



= s + s , 
A' = A(s', S12, S34) 



A = A(s:,-si2,-S34) = s'J + 2s'_(r + . (D.18 



Finally, a short remark on the integration of phase space in the above 2 — 4 and 2 — 5 
parametrizations is due: As all angular phase space variables are integrated analytically, 
and as the remaining invariant masses are integrated numerically without the use of Monte 
Carlo techniques, high precision numerical results can be obtained. 
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Appendix E 

The Cross-Section: Matrix 
Elements, Details, and Techniques 

In this appendix, matrix elements for process ( |1.1D will be presented. This includes the 
Born, the initial state bremsstrahlung, and the virtual initial state QED matrix elements. 
Together with the matrix elements, the analytical integrations of the angular variables 
are outlined. Results for the differential cross-sections obtained from these integrations 
are presented. Unless stated otherwise the ultrarelativistic approximation is used, that is 
the electron mass is neglected wherever possible. Correspondingly, final state masses are 
neglected too. A short account is given of the soft photon resummation technique which 
was applied to the 0{a) results. Generalization of the below results to process ( |4.2| ) means 
to add photonic Feynman diagrams with equal topology, as explained in chapter |. The 
necessary changes to couplings and propagators are most efficiently implemented with 
equation ( |4.^ ). 



E.l Born Level Results 

At Born level, the two Feynman diagrams given in figure |2.1| represent the amplitude for 
the pair production process ( |1.1| ). According to the Feynman rules from appendix [B.2| , 
using the four-particle phase space parametrization as in appendix |D.1| , and neglecting 
the masses of the initial state electrons and final state fermions, the Born t-channel matrix 
element is 



M 



B 



with 



a(3 



(E.l) 



u{pi) ■ [Lfjz-{1 + 75) + - 75)J ■ m(-Pj) 

^(-^2) ■7"[Le-e+Z-(l+75) +^e-e+Z-(l-75)] 'it X 

2ui-h) ■ 7" k-e+z'-(l + 75) + i?e-e+z'-(l - 75)1 " 7^ " uih] 



(E.2) 



62 



Spin indices are suppressed. The electron couplings to the Z° boson are 

Le-e+Z = , . ^ " ^ Sin26'^) 

4 sm&w costly ^ ^ 
^e-e+z = o • a ^ sin26'^ , (E.3) 

fermion couplings are derived from equation (|B.28| ), and the t-channel electron momentum 
is given by 

qt = ki-vi = -k2 + f 2 = ^iki-k2-vi + V2) 

= ki-pi-p2 = -k2 + P3 + P4 = ]j^{.ki - k2 - Pi - P2 + + Pi) ■ (E.4) 

Using the representation of four-vectors in terms of phase space parameters from ap- 
pendix |D.1| and the ultrarelativistic approximation (URA) one arrives at 

t = ql + ml ]^[A + Bcosd] , 

A = s - S12 - S34 , 

5 = . (E.5) 

The center of mass energy squared is given by s, the invariant masses are S12 and S34, 
■(9 is the scattering angle, and A = A(s, S12, S34) = s"^ + + -534 — 2ssi2 — 2SS34 — 
2S12S34. See appendix p.l| for further details. Finite boson widths are introduced into the 
Z^ propagators Dz{s) from equation 



Dz{si,) = Si^-Ml + i^Tzisij) , (E.6) 
where the Z° width is given by 

^'^''^^ = ^ ^ + ^^^') '^'^^^ ■ ^^-^^ 

with the color factor Nc{f). For the Born u-channel amplitude one obtains 

_ 9aa' gpp' ^ /M-/3' /M""' 

Bf = 2u{-k2)-r[Le^e+Z^-{l+l5)+Re'e+z'-{l-l5)]-<iul^ -Uiki) , 

Qu = ki-V2 = -k2 + vi = -{ki-k2 + vi- V2) 

= ki-p^-pi = -k2 + pi+P2 = ^{ki- k2+Pi+P2-p3-P4^) 
u = ql + ml '^=^ ^[A-Bcos^] . (E.8) 

It is appropriate to mention that the t-channel matrix element Alf and the u-channel 
matrix element Ai^ are obtained from each other by the interchanges 

vi^V2 , Mt2 ^ ^34 • (E.9) 



63 



After squaring the matrix element Ai^ = Aif , application of Feynman rule 8 from 
appendix p.2| , use of the phase space parametrization described in appendix p.l| , and 
fivefold analytical integration over the angular degrees of freedom of the four-particle 
phase space with the help of the integrals listed in appendices |G.1| and |G.2| yields the 
unpolarized cross-section 



J dsi2Pzisi2) J dSsiPzi-Ssi) X 



4mf 



Ami 



si2, S34) X 2 ■ BR{1) ■ BR{2) (E.IO) 

with branching ratios BR{1) and BR{2) of the Z° boson into the fermion-antifermion pairs 
/i/i and 72/2 respectively, and with 



Pz{Si 



vr \sij - M| + i,/sJ]Tz{sij)\'^ 



S{sij - Ml) 



A particularly convenient form for af{s; Su, S34) is obtained by using L^. 
= -(1 - 2sin^6'^) and Re = Re-e+z /{gzg2) = 2sin^6'^: 



^e-e+Z 



I i9zg2) 



[s; S12, 534) 



St^^i^] S12, S34) 



vr s 

^2 



g + (gi2 + ^34)' 
s - S12 - S34 

512,534) = 



-Cb-2 , 

1 Yn ^~ '^-^^ ~ -^34 + 
y/X s - S12 - S34 - 



(E.ll) 



^^^"(■s; S12, S34) decomposes into three components stemming from the t-channel the u- 
channel and their interference: 



Gl + g: + GT 



Gl 



Gi 



g: 



8S12S34 
1 

8 S12S34 



^ + 2s(Si2 + S34) - 8S12S34 + 4Si2S34(s-Si2-S34)/:b 
D 

^ + 4s(si2 + S34) + 16 s S12S34 [1 I Cb 

6 \ S — S12 — S34/ 



(E.12) 

To obtain the analytical result in equation ([E.ll| ), use was made of the representation 
of the processes Lorentz invariants in terms of phase space variables as outlined in ap- 
pendix p.l| . Throughout this thesis, algebraic manipulations were carried out with the 
help of SCHOONSCHIP [|0l, FORM [lOl, and Mathematica [|109|] . 
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E.2 Virtual Initial State Corrections 



Applying on-shell renormalization, the Feynman diagrams for the virtual initial state 
radiation to process ( |1 . 1| ) are shown in figure \i.2\ The "amputated" diagrams presented 
there must be connected to both the t- and the u-channel, resulting in a total number 
of eight diagrams. For our calculation, it is sufficient to use the t-channel virtual matrix 
element and compute its interference with the t- and the u-channel Born matrix elements. 
The interferences of the virtual u-channel with the Born matrix elements are then obtained 
by symmetry arguments as outlined at the end of this section. The virtual t-channel 
matrix element is obtained by substituting B^^/^qf + ml — ie) V^^ in equation (|E.1| ): 



M] 



vr^ 



(E.13) 



Using p as the loop photon momentum and applying dimensional regularization [|110| , |111 
one finds 



yap 



yap 

167r2 



- 1 fl 



(4-n) 



aP 



t,box 



(E.14) 



The counterterm part V^^^ given by 



v; 



ap 



t,CT 



b: 



aP 



qt + ml ■ 

aP 



— ( C^i + Cy2 + c. 

le ^ 



self 



qt + ml 



— f 2P + 4P^^ - 4 



(E.15) 



The counterterms C^i and 0^2 for the vertex corrections and Cseif for the electron prop- 
agator self energy have been derived in appendix |B.3 . 



Cself 



Cv2 



167r2 



6Zj 



2P + 4P^^ - 4 
-2P - 4P^^ + 4 



(E.16) 



The ultraviolet and infrared poles P and P are given in equations ( p.7| ) and ( p.8| ) of 



appendix p.5| . The matrix elements V^^'i, Vt°'y2^ ^"^eif^ V"^^^ belong to the positron 



aP 



t,v2i ^t,selfi 



,aP 



and the electron vertex corrections, to the intermediate electron self energy, and to the 



box correction as depicted in figure 3.2, In the ultrarelativistic approximation one finds 



aP 



t,vl 



^raP 



u{-k2 



pz 



7m it 7 



.p 



(p - qtY + ml-ie {p^ + 2/c2p) 



qf + ml — ie 



7 <it r 



(2A;i,^-^7m) 



qtY + ml 



le 



{jp' — 2kip) 



qt + ml 



^e-e+Z^-(l +75) +^e-e+Z^-(l -75) ' u{k 
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Taking into account that the virtual u-channel matrix element M.]^ is obtained from 
AiY with the help of the symmetry relation ([E.9|) one can now write down the virtual 
0{a) corrections to process (|1 . 1|) as the interference between the Born and the virtual 
matrix elements. Using the four-particle phase space parametrization from equation ( p.l|) , 



appendix p.l| , one arrives at 
(27r)^ 



an.) 



2s 



spins 



{Mr + -Mr) {Mf* + M 



B* 
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The evaluation of equation (|E.17|) is carried out in four steps. First step: The final state 
fermion decay angles are integrated with the invariant tensor integration formula ( |G.1| ) 
from appendix |G.1| . Second step: The result of the first step is used as input for a computer 
algebra program written in FORM ||108||. Inside the FORM program, steps three and four 



are carried out together with measures to simplify complicated intermediate expressions. 
Among the tools used to achieve simplifications, partial fraction decompositions played 
an important role. Third step: The photon loop momentum p is integrated over, making 
use of the integrals listed in appendix |G.5| . In the result of the third step, ultraviolet 
divergences from loop integrals and counterterms are explicitly cancelled. Fourth step: 
The integration over the boson scattering angle is performed with the help of the 



integrals given in appendix |G.6| . As a result of the above four steps, one can express the 
cross-section by a convolution similar to equation ([E.10|) 
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with a = e^/(47r) and 
lo = ln(s/mg) 



ivr 



ivr 



(E.18) 
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^virtual is identical to the 0{a) radiator obtained for the virtual initial state QED cor- 
rections in fermion pair production |5^, S^^^^^i contains all leading logarithms of 
the virtual cross-section. The infrared pole in S^^^^^i is cancelled by the soft part of the 
0{a) bremsstrahlung contribution. The 0{a) non-universal virtual cross-section does 
not contain leading logarithms I/3, but is rather complicated. It contains four parts corre- 
sponding to the four different matrix element products in equation ( PH.ITI) , 
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- 9s 
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In equation ( pi).2CI| ) the following notations are used: 

S — S34 



di2 

hi± 
h2± 



^B2 



S12 
s - S12 

1 



(compare appendix G.6) 
(compare appendix |G.6|) 
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Td- = L13 - L13 + L13 - L13 



S - S34 / V S - S34 / V S - S12 / V S - S12 

tmnrr. \ ^ . / tmin 



Tt. = Li3 -— ^ - Li3 - 



+ Li3 (- ;t%^) - Li3 (- ) (E.21) 



See appendix |G.6| for the definitions of Cb-, -^12, and £34. lyiq is found in integral 22.) of 
appendix p.5| . It is noteworthy that the quantities £b2 and Lbz have zero hmit for A ap- 
proaching zero. This is easily verified by expanding £5 in powers ofv/A/(s-a). It is seen 
from equation ( |E.2CI| ), the cross-section parts stemming from identical channels are much 



simpler than those originating from different channels. More specifically, o^nonuni ^^'^ 
(^I^nLuni do not contain trilogarithms whereas a^^nonuni and a^^nLuni do. 

A word about the relations between the above twofold differential cross-sections is in 
order. From inspection of the virtual Feynman diagrams one finds 

V,tt / \ V,uu / \ 

'^4,nonMml'^5 ■^12) -534^ — ^ i,nonuni\^ ^ ^12) 

'^l^nonunil'S) -312, S34) = Cr][i^*„^„-(s, S34, S12) . (E.22) 

The equalities in equation ([E.20|) are due to the invariance of (rX'^onuni and <ji'^onuni under 
the interchange S12 ^ S34. This is not surprising, because final state masses are neglected 
in the matrix elements and (J^'nonuni ^^d cr]fi^"„^^^ are kinematical functions, resulting in 
kinematically equivalent bosons. When adding the contributions given in equation ( |E.2CI| ) 
to obtain (lYnonuni 5 some minor simplifications of the overall expression are possible, but 
the representation ( |E.2CI| ) was chosen for reasons of transparency. 
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E.3 Initial State Bremsstrahlung 



The amputated Feynman diagrams for initial state bremsstrahlung are presented in fig- 
ure |3.1| . They must be connected to both t-channel and u-channel. Once more, the 
symmetry of the problem simplifies the calculational task. It is sufficient to compute the 
square of the t-channel bremsstrahlung matrix element and the bremsstrahlung t-u inter- 
ference. The u-channel contribution to the cross-section is obtained from the t-channel 
contribution by virtue of a symmetry operation similar to relation ( [E.22|) . The brems- 
strahlung t-channel matrix element is derived from equation ( |tL).l| ) via the substitution 



Bt I i^t + "^e — i?" £^ with polarization index A. Neglecting numerator electron 

masses, the Feynman rules in appendix [B.2| yield 
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(E.23) 
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(E.24) 



The momentum of the radiated photon is called p. Similarly, one finds for the u-channel 
bremsstrahlung matrix element 
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Z2 Ul J 
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2 3 

— T 
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In the above matrix element the radiatively changed Mandelstam variables were intro- 
duced. In terms of the five-particle phase space parametrization from appendix |D.2| they 
are given by 



h = 
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{ki 
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- hi cos9f( — c sm9RC0S(f)R 
a2 — 62 cos9r — c sin6'i{ coscpR 
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Z2 



-2 hp = y(l + /5cose) = 
-2k2P = y(l-/3cos^ 



■Zl 



(E.26) 



with the kinematical parameters 
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derived from the representations ( |D.10| ) and ( p.l5| ) in appendix |D.2| . The Lorentz invari- 
ants satisfy the relations 

ti + Ui + Zl = s — a 

t2+U2 + Z2 = s - a 



Zl+ Z2 = S_ 

The total bremsstrahlung cross-section is obtained by the integration 



2s 



(E.28) 



(E.29) 



spins 



with the five-particle phase space dV^ as parametrized in equation (p.9|) , appendix |D.2 . 



After invariant tensor integration over the boson decay degrees of freedom with for- 
mula ( |G.1| ) from appendix a FORM [p.08|| computer algebra program is used to 



further evaluate equation ([E.29|) . Inside the FORM code, partial fraction decomposition 
is used to simplify intermediate expressions, and the following non-trivial scalar products 



are needed: 
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The, after the tensor integration step, remaining three angles Or, (pR, and finally 6 (see 



appendix |D.2| ) are integrated analytically with the help of the integrals listed in appen- 
dices |G.3| and |G.4| . The complexity of this threefold analytical integration is due to the 
appearance of the angularly dependent Mandelstam variables ti, t2, ui, U2, Zi, and Z2 in 



the denominators of the matrix elements TWf'^ and M^'^, equations ( [li].23 ) and ( [ti].25| ). 



The result for the total bremsstrahlung cross-section splits into a universal, leading loga- 
rithm part with the Born cross-section factorizing and a non-universal part. 
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(E.31) 



The non-universal cross-section term cf^nonunj is a sum of the t-channel, the u-channel 
and the t-channel/u-channel interference bremsstrahlung contributions. 
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R,tu 
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The shorthand notations L^i, Lc2, Lcs, Lc4, Lc^, L^q, Lcj, Lcs, D\,Dl,, DU, -D^^i, D 

Dz2tl, Dz2t2 , -D*", D% Dtiu2, Dt2ul, Dl^u2-> ^t2ul^ A'l«2 5 Df^^ 

pendix |G.4| . For easier reference, the number of the integral from appendix |G.4 



125 ^zltl, J^zU2, 

are introduced in ap- 
corre- 
^ from 



sponding to each of these notations is given in table Together with a 
equation ( |E.2CI| ), cr^nonuni from equation (Pi).32|) represents a main result of this thesis. 



Notation 


Lci 


Lc2 


Lc3 


Lc4 


Lc5 


-^c6 


Lc7 


Lc8 


Number from 
appendix |Cj.4| 


10.) 


13.) 


30.) 


31.) 


10.) 


26.) 


27.) 


26.) 


Notation 


D\ 


Di 
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-D^iti 


Dzlt2 


Dz2tl 


Dz2t2 




Number from 
appendix G.4 


18.) 


19.) 


46.) 


10.) 


11.) 


12.) 


13.) 


26.) 


Notation 




Dtlu2 


Dt2ul 


^tlu2 


^t2ul 


^tlu2 






Number from 
appendix p.4| 


27.) 


51.) 


52.) 


53.) 


54.) 


57.) 


58.) 





Table E.l: Shorthand notations from equation ( \E.3^ ) together with the number of the 
integral from appendix \G.^ where each notation is defined. 
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The radiator (s^ + s'^) / {s'^s'_) obtained in equation ( [til.31| ) is known from e~^e~ annihilation 
into fermion pairs ||112|| . It is seen from equation (p.ll|) that the radiator's singularity in 
the factor = l/(s — s') oo is related to the radiation of soft photons. This 

so-called infrared singularity is identically cancelled by the infrared singularity arising 
from the virtual corrections. The cancellation of infrared singularities is worked out in 
appendix ^A . 

It is important that the non-universal contribution cr^onuni is not infrared divergent. This 
was explicitly checked by taking the hmit s' ^ s for all parts of (T^nonuni that contain 
factors \/s'_. In addition it was checked that all parts of o"f„o„„„j have finite limits for 

^0. As a concluding remark it is mentioned that, similar to the case of virtual 
non-universal corrections, some simplifications are possible when terms are collected from 
the individual contributions in equation ( |E.32| ) to yield cr^nonmi- these are not major, 
the expression for crfrionuni not rewritten. 



E.4 Treatment of Infrared Singularities 



To show how the infrared divergences from the virtual and the real bremsstrahlung cor- 
rections cancel, a short derivation of this cancellation is worked out below. One starts 
from the bremsstrahlung matrix elements given in equations ( |E.23|) to ( [E.25|) . The in- 
frared divergence originates from the square of those parts of the matrix elements which 
contain the bremsstrahlung photon momentum p in the denominator. In equations ( |E.23|) 
to ( [ti].25| ) these are the parts containing 1/2:1 and 1/2:2 as is seen from equation ( [ti].26D . 
Denoting the matrix element parts yielding the divergence by using very small p , 

i.e. neglecting numerator terms proportional to p, and applying the Dirac equation one 
gets 
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IR 



ki.e k2.e 
ki.p k2.p 



(E.33) 



with the Born matrix element M.^ = M.f+Ai^ taken from equations ( [E.ID and ( [E.8| ), and 
with the photon polarization vector e. Evaluation of the fully differential bremsstrahlung 
cross-section from and summation over photon polarizations yields the factorization 

(E.34) 
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(E.35) 
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with the fully differential Born cross-section da^ , the photon energy p'', the soft cutoff 



A>0, the Heavyside step function 0(a;), and the long known Low factor ||113 



■IR 



ki k2 
ki.p k2.p 



(E.36) 



The infrared divergence is obviously contained in ^^o/f Following a procedure of dimen- 



sional regularization similar to the one outlined in reference ||llCm and carried through 
in 



114 one arrives at 
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A derivation of S^^jt is given at the end of this section. The result may be verified against 



reference |63] where regularization with a photon mass A is used by the replacement 
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The dependence of 5^^/* '^^ the cutoff u is cancelled by a similar dependence in Sj^^.^. 
Integration of equation ( [E.34[ ) over the Born phase space yields the infrared divergent 
cross-section contribution with the factorized Born cross-section 
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Next, the contributions from 0{a) virtual and bremsstrahlung corrections must be added 
to cancel the infrared divergence P^^. As only universal cross-section parts are infrared 
divergent, it is sufficient to consider those. Non- universal contributions are therefore 
ignored for the remainder of this section. First, the universal cross-section contributions 



from equations ( pi]. 18 ) and ( E.31| ) are rewritten using the definition L = ln(s/ 
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One observes that addition or subtraction of non-divergent terms in the second integrand 
of equation ( [E.41| ) yields neghgible effects for sufficiently small infrared cutoff a;. Therefore 
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It is evident that the infrared divergence of the bremsstrahlung part is in the second term 
of equation ( [li).42D . The infrared divergence is located at s = s' and the corresponding 
cross-section is given by 



cr 



IR 

soft 



a f ds' 2 s B , 

-■ j-L-a [s] 

TT J S S — S' 



(E.43) 



Replacing this expression by the expression given in equation ( [E.39| ), one can rewrite 
equation ( |E.42| ): 
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The complete universal 0{a) cross-section is obtained by adding the virtual and the 
bremsstrahlung parts from equations ( |E.4CI| ) and ( Pj^.44| ). The result, also known from the 
literature on fermion pair production in e~^e~ annihilation p.l2|] , reads 
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with S^^f^^i from equation (|E.18| ) and 6^^j^ from equation (|E.371 ). The dependence on the 
infrared cutoff u cancels if the integration is carried out. From the representation ( |E.45| ) 
of the universal 0{a) cross-section the cancellation of infrared divergences is obvious. 
As stated earlier, the non-universal 0{a) cross-section parts are infrared finite. One 
can therefore present a singularity-free full 0{a) cross-section by adding non-universal 
0{a) virtual and bremsstrahlung contributions from equations ( [E.18 ) and (|E.31|) to the 
result (F^ )- 
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Derivation of the soft bremsstrahlung infrared divergence 6^^. 



soft 



To derive the expression for 5^^^^ given in equation ( [E.37| ) on starts from equation ( [E.35| ). 
Introduction of a Feynman parameter a yields 
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Using the spatial angles 9i, 
ki, the above denominators are expressed as squares of 

— ki-p = k^p^ (1 — /?j cos6'j) with velocities 
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1, 2, a between the photon momentum p and the momenta 



A = 1^ • (E.47) 



To dimensionally regularize equation ( [li).46D , the transition from 3-dimensional to {n — 1)- 
dimensional space has to be performed. As all three terms in the integrand of equa- 
tion (|E.46| ) do not depend on the polar angle 9 of the photon emission, one has the 
freedom to choose a coordinate system with 9 = 9i for each of the three integrand terms in 
eq. (|E.46| ). Physically this freedom of choice is equivalent to the isotropy of soft photon 
emission. With the mass scale /x to maintain the dimension of the expression, one gets 
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Introducing the integrals 
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sm 



^''-^'^9,d9. 



n — 4 n 



+ 1 n — 4: 

'1 — x^) 2 dx 



1 + (n-4) In— + C(n-4) 
/i 



(3i cos9i 



(1 - P^X) 

dx 

[ {1-P.xf 



-1 
+1 



1 + 



n 



ln(l -x^) + 0(n-4) 



(E.49) 



neglecting terms of 0{n—A), and taking into account that /3i = /32 = /5 = y 1 —Ami/ ^ 5 
finds 

1+1 ,1+1 



eIR 

"soft 



+ In: 



A 



x) 



■ Jda Jdx J-'{a,x) + Jdo^ Jdx \ia{l — x^) J-'{a,x) 

-1 

s — 2ml 





Ami 



Ami 



(A;0)^(l-/?„x)' s{l-Pxy s{l-Pxy 



(E.50) 
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Using the two integrals 

7 1 

'dx 

-1 
+1 
'dx 



-1 



{l-^3^xf 

ln(l - x2) 



2 



21n2-lln-^+^^ 



A 1 - A 



(E.51) 



results in 



1 +1 



2 jdoi Jdx J-'{a,x) 

-1 

-Jdajdx ln(l-x^)^(a,x) 



i 



da 



- 2 



ln2 



s — 2m? 



da 



- 2 



+ 



da , 1— /3q 
m- 



1 1-/5 
(E.52) 



with the Feynman parameter integrals 



da 

_h.2 



da 



[s — AmD a (1 — a) + m 



A 1 

s/3 ^1-/3 



m- 



■ In- 



i-N{i-2ay 



■In 



1-Px 



s(3 J x{l-(3x){l+(3x) l+(3x 



2 

Jp 



2Li2(-/3) - 2Li2(/?) + Liaf 



Ll2 



2 2 



-In- 



2 y 

(E.53) 

■1-/3 -,nl - • t^-^^' 

taking the result to the ultrarelativistic limit, and evaluating it at the scale fi = ml yields 
exactly the result given in equation ( [E.37| ), 



2 2 

Introducing the results (|E.52| ) and ([E.53|) into equation ( [E.50| ), using 



1 + /3 URA , S 

In = m — ; 



"soft 









2- 











7r 



+ Inu; ( In-- - 1 1 + ^ In^-- - — . 



TT 



mi 



(E.55) 
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E.5 Soft Photon Exponentiation 



The factor iii the universal part (see equation ( [E.45| )) yields large contributions at 
the phase space edge close to s = s' . Therefore resummation of higher order contributions 
is suggested to obtain reasonable precision. One can follow the derivation of the resummed 



cross-section given in reference ||6^ , because its authors treat precisely the coUinear leading 
logarithms that are found in the universal cross-section parts of the presented calculation. 
To connect to the formulae of ref. |63|, one starts from equation ( |E.45|) and introduces 
the notations 



V = 1 



Si = 
Hi = 



2a 

7r 
a 

TT 

a 

IT 



S 
S 

■ L 



5. 



2\nuj-L 



Heis'/S) 



a 

7T 



7!" 
"3 

1 r 



2-v 



(E.56) 



From equation ( [E.45| ) one obtains the cross-section for process (|1 . 1|) with universal initial 
state QED corrections 



a 



ISR 



dv a^is^ 




3^ 

1 + Pe^nu + Si + ^ In^u; + /?e Incu ^1 + ^2 



Pe - + Hi + P'e— + Pe- Si + H^ 
V V V 



+ 0{a' 



(E.57) 



with the second order initial state radiation taken from reference |^3| where expressions 
for the second order virtual+soft part 5*2 and the second order hard part H2 can be found. 
Before carrying on, it is worth mentioning that the infrared divergent second order soft 
initial state radiation part S^^j^ ,.^jr^ is given by 



^IR 



2 "soft 



+ C2 



(E.58) 



with infrared finite C2. This remarkable property indicates the possibility to obtain higher 
order initial state soft-coUinear divergences from an exponentiation of the lowest-order 
apart from a finite correction factor c„ K]. 



correction 



Equation (Pi).57|) indicates that the cross-section will be enhanced by terms ln"~^f /v 
which become especially important near the phase space boundary v ^ 0, i.e. the en- 
hancement is due to soft photon emission. It will be present in all orders n of perturbation 
theory, and therefore an appropriate resummation technique is needed. As the logarithms 



originate from infrared singularities they can be resummed |63], and, according to the 
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Bloch-Nordsieck theorem ||115|| and Yennie, Frautschi, and Suura ||116|| , they indeed ex- 



ponentiate to all orders of perturbation theory. Therefore equation ([E.57|) represents the 
truncated power series expansion of the universal cross-section to process ( p. . 1| ) with soft 
photon exponentiation which reads 



a, 



ISR 

exp 



J dv a^(s') [5{v) ■ S ■ e^=^°" + Q{v - uj) ■ h] (E.59) 



with the regular soft+virtual part 5* and the hard part if, 

oo 

S = '^Sn with 5*0 = 1 for the Born contribution , 

/5e" In"-^ V 



n=0 



H 



E 

ra=l 



{n-l)\ 



■S + Hr, 



(E.60) 



All Hn are integrable in v, and the infrared divergent term in H is completely determined 
by the requirement that the Incj terms appearing in the v integration of H cancel all 
Inu; terms in the virtual+soft part S-e^'^'^^. To rid equation (E3^) of the unphysical 
parameter oj one introduces H = J2'^=i Hn , uses the equality 



E 



In"-^ V 



n=l 

and rewrites 



n-l)\ 



l3e-l 



(EM) 



a, 



ISR 

exp 



dv (T^(s') { \6{v)-e'^^^'''' + Q{v-uj)-(3ev'^^-^]-S + Q{v~uj)-h\ . 



For the three above integral terms one finds 



(E.62) 



(E.63) 



dv a^{s')-Q{v-uj)-f3eV^^-^-S 



J dv a^{s')-PeV^''^-S - j dv (r^{s')-PeV^^-^-S 



and, because H is integrable and u is sufficiently small, 



j dv a''{s')-e{v-uj)-H = J dv a''{s')-H . 





(E.64) 



(E.65) 
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Adding the results of equations ( |E.63| ) to ( |E.65| ) one notices that the contributions 
from ( [E.63|) and from the second term of (|E.64|) cancel each other. Thus a new and 
physical form of the initial state radiation corrected cross-section with soft photon expo- 
nentiation is obtained from equation ( p.62| ): 



a 



ISR 

exp 



J s 



(E.66) 



A cross-section with soft photon exponentiation is now obtained by replacing the 

terms S and H in ( [E.66D with corresponding truncated sums from equation ( [E.60D . For 
an 0{a) ISR corrected cross-section with soft photon exponentiation the series' (|E.6CI|) 
are truncated after the first order term which means 



S 
H 



Hi + Oia^) 



Returning to process ( |1 . 1| ) and aiming at the 0{a) ISR cross-section with soft photon 
exponentiation, it is necessary to include non-universal cross-section parts. Since the 
virtual and the bremsstrahlung non-universal parts as derived in appendices pi).2| and [E.3 



are both infrared finite and do not contain any singularity, they may be readily added to 
the regular first order expressions Si and Hi without changing the argument leading to 
equation ( [E.66|) . Therefore accounting for 0{a) non-universal cross-section parts means 
to replace: 



1 + Si 



a''(s')-il + Si] + a, 



ct''(s')-Hi 



<y''{s)-Hi 



+ s 



V 

nonuni 



ds' 



Si 
Hi 



(E.67) 



with CTnanuni fro^i cquatiou ( |E.18| ) and CTnonuni taken from equation ( [E.31| ). Thus one 



obtains for the 0(a) ISR corrected cross-section with soft photon exponentiation for 
process 



a. 



ISR,0{a) 
exp 



ds^ 
s 



-'■Si + Hi} . 



(E.68) 
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Appendix F 

Evaluation of Loop Integrals 



Calculations of radiative corrections within the Electroweak Standard Model contain loop 
integrals as basic building blocks. For the one-loop case, loop integrals have been inten- 
sively studied by several authors |104| , |117| , |118|| , and the literature also contains 
pedagogical treatments |9^. It is noteworthy that scalar one-loop diagrams with five 
or more external lines can be reduced to the scalar 1-, 2-, 3-, and 4-point functions [p.04 



As tensor m-point functions are reducible to linear combinations of scalar fc-point func- 
tions with k < m, it is sufficient to know all possible scalar 1-, 2-, 3-, and 4-point functions 
together with the appropriate tensor reduction rules. As it is not the intention of this 
appendix to repeat the results obtained and comprehensively presented earlier (see e.g. 
references p6| , |118| ]), it will be restricted to the formulae relevant for this thesis. This 
means a restriction to the 1-, 2-, 3-point functions, because the only 4-point function 
that had to be dealt with could be easily evaluated by direct integration. It should be 
noted that, in special cases, the more general results obtained below may be subject to 
simplifications by virtue of the relations presented in appendix |C|. Such simplifications 
can often be more easily seen, if the specific loop integral is directly evaluated. In this 
case, however, the general result is still very useful for cross-checks. Subsequently, the 
following quantities will be investigated: 



Aoiml) ; 



-(n-4) f _}_ 



-(n-4) 



(27r)" Di 

{1, g^, q^qy] 



(27r)" D^D2 
d'^Q {1, qf,, q^^qu, qf^quqp} 



(27r)" 



Di D2 Ds 



(F.l) 



Dimensional regularization is applied to render all expressions well-defined. is an arbi- 
trary mass scale. The following definitions are used: 



q^ + ml — is 



(q + Pi) 



' + m\-ie 



(F.2) 



The n-dimensional Dirac algebra needed for the evaluation of loop expressions can be 



found in [^, |119|| or in standard textbooks such as [0. Care is required if 75 matrices 
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are involved. Note that the expressions in equation ( [t*'.!! ) are invariant under shifts of the 
loop momentum q. Before proceeding further a lemma will be presented that is going to 
be very useful in section p^'.l| . 

Lemma 1 Let a he real. The numbers b,c and yo may be complex. Let further —e be 
an infinitesimal real number carrying the sign of I = '^m{ay'^+by+c) , where I is supposed 
to not change sign for ye [0,1]. Let the infinitesimal real number —6 carry the sign of 
'^m{ayQ + byo + c) . Let yi and y2 be the roots of ay'^ + by + c = 0. Then the identity 
1 



dy 



y-yo 



In 



ln(ay^ + by + c)- \n{ayl + byo + c) 
?/o - 1 



v{-yii-y2) - v{yo-yuyo-y2) - 77(0 - ie, a + 15) 



+ R(2/o,2/i) + ^(2/0,2/2) 
holds, if the two parameter function R is defined by 



R(2/o, yi 



dy 



y-yo 



yo 



ln(2/ - yi) - ln(yo - Vi) 
1 \ 



,yo - Vi 



Li2 



Vo 



+ v\ -yi, 



Vo - yi, 



■In- 



.2/0 - yi 
yo 



yo - Vi 



- ^ l-2/i, 



1 



■In 



2/0 - yi, 



Vo - 1 
2/0 - Vi 



(F.3) 



(F.4) 



The proof of Lemma |l| can be found in appendix B of reference ||104|| . See appendix |C.l 
for the definition of ri{a, b). 



F.l Scalar m- Point Functions 



The evaluation of scalar m-point functions in n dimensions is based on the represen- 
tation of n-dimensional space in spherical coordinates as outlined in the appendix of 



reference ||110|| . For m > 2, the second important ingredient is Feynman parametrization. 
One arrives at expressions which become rather handy for n = 4 — e, |e|^L Integrals 
converge for e > in the case of an ultraviolet (UV) divergence and for e < in the case 
of infrared divergences. 

The scalar 1-point function is given by 



Ao{ml] 



1 



A 



16 7r2 
2 



m 



A - 1 + In 
+ 7b - ln(47r) , 

n — 4 

with Euler's constant '-/e = 0.577216... . 



m| — le 



(F.5) 
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The scalar 2-point function is given by ||104 



16 7r2 



A + In' 



-Pi - ig 

/i2 



+ E 



1=1 



ln(l — Xi) — Xi-hi 



Xi - 1 



- 1 



(F.6) 



where Xi and X2 are the roots of —pfx"^ + {pf + m| — m^)x + ml — ie = . If = 0, 
there is only one root, and ln(— — ie) should be replaced by ln(m2 — ml — ie). If in 
addition ml = m\ = m^ one finds 



5o(0, m^m^^ 



1 I ^ , — le 
A + In — 



16 TT 



(F.7) 



The pathological case p\ = m\ = m\ = is infrared divergent and would require appro- 
priate special treatment. 

The scalar 3-point function has a more complicated structure, but is not ultraviolet di- 
vergent. It can be presented as 



1 X 



CQ{pi,P2,ml,ml,ml) 



16 7r2 



16 71' 



dx / dy ax^ + by^ + cxy + dx + ey + f 







Co 



(F.8) 



with 



a = -pI, b = -pI, c = -2pip2 =pI+pI- {pi + P2)^ 

J 2 2,2 2 2|2io J! 2- 

d = m2 — m^ + P2, e = m^ — m2 + Pi + 2pip2, j = m.^ — le . 
Defining a as one of the roots of ba^ + ca + a = Q , the first integration yields: 



Co 



dy 



- dy 



In 



by'^ + {c + e)y + a + d + f - \Yibyl + {c + e)yi + a + d + f 



(c + 2ab)y + d + ea + 2a + ca 



(l-a) 



(c -|- 2ab){l — a)y + d + ea 



In 



[a + b+c)y'^ + {d+e)y + f 



In 



[a + b+c)yl + {d+e)y2 + f 



+ 



1 a ■ I In ay"^ + dy + f 
J (c + 2ab)c 



In 



ayl + dy3 + f 



+ 2ab)ay — {d + ea) 



(F.9) 



with 
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yi 

2/2 
2/3 



yo + a 

yo 
1 — a 

_ yo 

a 



yo 



d + ea 
c + 2ab 



Taking into account that for each term of equation ( |b'.9| ) the non-logarithmic factor can 
be written as (c + 2ab) ■ [y — yi] and using the identity ( |F.3|) presented in Lemma || at 
the end of the previous section, it is straightforward to further integrate equation ([F.9|) . 
This solution is valid for all physical cases, i.e. for physical four-momenta pi and p2 and 
physical masses ||104|| . Physical four-momenta are real, masses may be complex but must 
have a positive real part. 

For real roots a and stable particles, the 3-point function can be cast into a rather simple 
form: 



a, rrii e M 



3 2 



with 



Xi 



X2 = - 



1 y y (-1)'+' 



+ 2a + (e + c)a 
c + 2ab 

d + ea 



Li2 



Xi 



xi - yij 



Lis! 



Xi 



1 



(F.IO) 



x^ 



1 - a){c + 2ah) 
d + ea 



a{c + 2ab) 



yij 
y2j 

ysj 



c - e ± y'(c + e)2 - 46(a + d + f) 
2b 



-d-e±^{d + e)2 - 4/(a + b + c) 

2{a + b + c) 
-d ± y/d' - 4a/ 
2a ■ 



Among the special cases of the scalar 3-point function, one is rather interesting, namely 
the infrared divergent case which is encountered, if one propagator mass vanishes and the 
momenta pi and pi + p2 are on-shell. Without loss of generality one can choose 

{Pi+P2y = -ml 



pI 



-m 



2 ' 



A- 



mi = , 

The infrared divergent scalar 3-point function is then given by 
Co^(Pi,P2,0, m^,m^) = 

-(n-4) 



(F.ll) 



1 



(F.12) 



(27r)" (g2 - ie)(g2 + 2qpi - ie){q'^ + 2q{pi + P2) - ie) 

From power counting in equation (|F.12|) it is seen that is infrared divergent for 

n = 4, but regularized by n > 4. Expansion around n = 4 yields 

^IRy^ ^ n ™2 



Co (Pi>P2,0,m2,m3) 



1 



16 7r2 



dy 







-Pi 



A^^ + In 



:F.13) 



with —Py = —p2y'^ + {p2 — m2 + ml)y + m2 — ie and A^^ = A, the superscript IR indicating 

^2 



its infrared origin. The case P2 = yields a simplified integral and will not be treated 
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here. In addition, the cases 1712 = and = are excluded from further consideration 
to ensure that the roots yi and 1/2 of —py = do not coincide with the integration edges. 
Defining 



4 

2 



one can rewrite 



2 2 I / 2 2 I 2\ I 2 

2 2 I / 2 2 I 2\ I 2 

-P2Z/2 + {P2 -^2 + m3)y2 + m2-ie 



Cl^iPuP2,0,mlml) 



16 7r2 2 -Vr{yi-y2) 



dy 



y-yi 



A^^ + In ■ 



-P 



yi 



+ In(-p^) - ln(-p^ 



Using 



dy 



y-y2 



^IR + In ^ + _ ln(-p: 



-P 



(F.14) 



In 



2/i - 1 



(F.15) 



and identity ([F.3|) from Lemma |l] further integration is straightforward. Despite the 
generahty of this result direct integration of the R.H.s. of equation ( |F.13| ) is simpler in 
many practical cases. 



F.2 Reduction of Tensor Integrals 

Because of Lorentz covariance, tensor loop integrals can be decomposed into linear com- 
binations of all possible tensors with identical rank that can be built from the involved 



external four-vectors and the metric tensor. 

= Pi^B, (F.16) 

B^^ = Pit,PivB2i + gt,vB22 (F-17) 

C, = Pl,C^l + J92^Ci2 (F.18) 

C^u = Pl^PluC21 + P2iiP2uC22 + {piP2}f,uC23 + g^^uC24: (F.19) 
C^^^'P = PltJ,Pl^PlpC-ii + P2fMP2uP2pC32 + {p 2P1P1} pup C 33 

+ {piP2P2}pupC34. + {pig}pupC35 + fi'}/.^p C36 (F.20) 
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The linear factors Bn^j-j and Cij are linear combinations of scalar loop integrals and depend 
on the Lorentz invariants that can be formed from the momenta Pi of the involved external 
four- vectors and on the propagator masses squared mf. For the sake of simplicity, no 
arguments are exposed in equations ( [F.16| ) to ([F.20| ). If nothing else is given, arguments 
are always assumed to be as exhibited in equation ( [b'.l| ). The convenient shorthand 
prescriptions 



{p 

{Pi9]^lup 



= Plf,P2u + P2,,Plv 

— PifJ-Pji^Pip ~l~ PifJ-Piy Pjp ~l~ PjlJ-Pif Pip 

— Pip- 9up ~\~ Piu 9pp ~l~ Pip Qpu 



:f.2ii 



were taken from reference ||118|| , the other notations are rather standard (see e.g. |p6| ). 
Now the computation of the linear factors and Cij is described. To evaluate Si, both 
sides of ([F.16|) are multiplied by pi. This is followed by quadratic supplementation of the 
integrand's numerator on the L.H.S. Then division by p\ yields the well-known result 



2 Pi 



fiBo + Ao{ml) - Ao{ml 



[F.22) 



with fi = ml — m\ — p\. Similarly, multiplication of both sides of ( [bM?] ) with g^^^ and 
PiPi respectively yields a linear system of equations. 




;f.23) 



which has the solution 
1 



B 



22 



n — 1 



2 



mlBo + ^Ao{ml] 



B2I = -o 

Pi 



^f,B^ + ^Ao{ml)-B22 



(F.24) 



The above already reveals the algorithm for the reduction of tensor n-point functions. Mul- 
tiplication of both sides of the tensor integral's linear decomposition (see equations ( \F.1(\ ) 
to l \F. 2(X j ) by the integral's tensor coefficients and subsequent evaluation of the resulting 
integrals with the help of quadratic supplementation of the integrand yields a linear system 
of equations for the linear factors Bn^j) and Cij in terms of scalar m-point functions. This 
resulting system originally consists of k equations for the k linear factors to be evaluated. 
A more elegant way with reduced dimension of the system of equations to be evaluated 
is given in references [^, |118|| . Subsequently, we shall follow their approach and present 



the solutions for the remaining linear factors Cij. First define 



B 



{0,p,pv) 



-(n-4) 



d^'q {1, q^^, q^,qu} 



(2vr) 



D, D, 



(F.25) 



with Qi = and Q2 =Pi. The denominators for this definition are found in equation ( [l''.2D . 
The arguments {k,l), k <l will be correspondingly used for the linear factors Bi(^j^{k,l) 



which can be derived by virtue of equations (|F.16|) and ( |F.171 ). Using /2 = — rrig + 
(pi + P2y and considering 



pj 
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X 



Pi P1P2 
P1P2 pI 



an invertible matrix (i.e. assuming {P1P2Y ~ P1P2 7^ O); ^^e linear factors Cij for the 
R.H.s. of equation ( [b'.18D are given by 



'a 
c 











Ri 




R2 


= P'2C, 



1 / Ri 
R2 



(F.26) 



/iCo + i?o(l,3)-i?o(2,3) 
/2Co + So(l,2)-5o(l,3) 



The hnear factors on the R.H.s. of equation (|F.19|) are given by similar matrix rela- 
tions, namely 



C21 

,C'23, 



X 



-1 R3-Q 

Ra 



24 



C23 
C22 



( R5 



The UV divergent linear factor 



24 



n-2 



^ICo - 7: (fiCu + /2C12 - 5o(2, 3 



is obtained from the original 5x5 system of linear equations and 



R3 
R4 

R5 
Rp. 



P1P2 Cf,u - pI R5 

P1P2 ^ 

{P1P2) ■ P2P2 ~ pI ■ P1P2 ^ 

{piP2f - P1P2 

{P1P2) ■ PiPi - pI ■ P1P2 



{pmf - P1P2 



P1P2 



/iCn + 5i(l,3) + 5o(2,3) 
/2Cn + 5i(l,2)-5i(l,3) 
/iCi2 + 5i(l,3)-fii(2,3) 

/2Ci2-fil(l,3) 



P1P2 Cf,u - pI Ra 

Finally, the matrix relations for the linear factors C^j are given by 

-1 / Rw \ 



C35 
Cm, 

C'si 



X- 



X- 



\Rii J 
1 / -R12 — 2 C35 



R 



13 



C34 
C32, 



i?14 

i?l5 — 2 C36 



In addition there is the redundant matrix relation 

— 1 / -^16 ~ C*36 \ 



C33 
,^'34 , 



X- 



Rl7 ~ 



(F.27) 



(F.28) 



(F.29) 
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The scalar factors needed are 



Rio = 


1 

2 


/iC'24 


+ -B22 




- ^22(2,3) 


Rii = 


1 

2 




+ -B22 




- 522(1,3) 


Rl2 = 


1 

2 


/lC'21 






-So(2,3) 


Rl3 = 


1 

2 


/2C2I 


+ ^21 


.1,2) 


- ^21(1, 3) 


i?l4 = 


1 

2 


/lC'22 


+ ^21 




-S2i(2,3) 


Rib = 


1 
2 


/2C'22 


— -B2I 


:i,3) 




Rl6 = 


1 
2 


/lC'23 


+ ^21 


;i,3) 


+ 5i(2,3)" 


Rl7 


1 
2 


72^23 


— -B2I 


:i,3) 





This completes our outline of the reduction of tensor to scalar loop integrals and the 
appendix on loop integrals as a whole. 
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Appendix G 

Table of Analytical Integrals 



In this appendix the analytical solutions of the integrals needed for the semi-analytical 
cross-section calculation of process ( |1.1| ) are presented. These integrals are necessary to 
obtain the twofold differential cross-section d'^a/dsi2ds34 in the Born and virtual ISR 
cases and the threefold differential cross-section d^^a/ds'dsuds^^ for initial state brems- 
strahlung. This list is restricted to integrals needed for the t-channel contribution and 
the t-channel/u-channel interference, because the contributions from the u-channel and, 
in the virtual ISR case, from the u-channel/t-channel interference are obtained by sym- 
metry arguments as detailed in appendix |^. Integration of the azimuth around the beam 
axis is trivial. The integrals presented below have been numerically checked. 



G.l Fermion Decay Angle Tensor Integrals 



Integrations over final state fermion decay angles are carried out with the help of invariant 
tensor integration. Using the decay matrix element M/^ from equation (pi).2|) one obtains 



E 



Spins i,j 



8 Si 



E 

Spins i,j 



dcos6. 



2iv 



8 s,: 



/VT / y 2 ^ 



with rnyi 



mi and 77134 



(yfj.jg and g^j^ are found in appendix |B.2 



?7i2. The vector and axial vector boson-fermion couplings 
The result of the above decay phase space 



integral factorizes in the cross- sect ion. Therefore one can use the above formula ( |G.1| ) 
to integrate the final state fermion decay angles for both the 2—^4 and the initial state 
bremsstrahlung 2^5 particle phase space. The factor [g^^j^ + dffB^^ enters the decay 
width Tz{sij), equation (|2.3|), and thus the Breit-Wigner density pz{sij), equation (|2.2[). 
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G.2 Born Level Phase Space Integrals 



For the calculation of the Born level cross-section of process ( |1 . 1|) only a few integrals over 
the boson scattering angle are needed in addition to those over the final state fermion 
decay angles. They are listed below, making use of the phase space parametrization given 
in appendix p.l| , the ultrarelativistic approximation (URA), and the notation 



+1 



d cost!} A 



(G.2) 



-1 



10 

2.1 



6.1 



[1]* = 
[cosi9], 



t2 



1 i^"5-^+V^ 
a/a S — cr— a/a 

1 



3.) [cos^^]^ 
5.) 



1/3 



In- 



-cr 



-a- 



Va 



G.3 Bremsstrahlung Integrals — First Series 



For the calculation of the bremsstrahlung process ( |1 . 1|) , the phase space was parametrized 
as exposed in equation ( |D.9[ ), appendix p.2| . After integration over the final state fermion 
decay angles, the first series of bremsstrahlung integrals is over the scattering azimuth 
and polar angles (pn and of the boson three-vector vi in the two-boson rest frame. 
Using the double integral notation 



^ 2n +1 



(G.3) 



and the ultrarelativistic approximation, i.e. neglecting the electron mass wherever possi- 
ble, the bremsstrahlung integrals for the first series are listed below: 



1-) [1]. 



2.) [cos0R ■ /(cos6'ij)] 



R 







3.) [cos6'h] 



R 



4.) [cos20^] 



R 



1/2 



5.) [cos^^h] 



R 



1/3 



6.) [cos^^ij, ■ cos'^Or 



R 



1/6 
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■1" 



R 



2s' 



A'(s^ — s'_ cos6 
2s' 



X' S, 



■I 



ti 



dl 



\n( ''+('''+V^)-2s'o-+si5-(s'+'5+A/V)sicos0' 



- 1 ■ 

T2. 



R 



2s' 



VX'{s[_ + si cos6 
2s' 



X' S, 



■I 



t2 



d2 



■In 



sUs' + ^/x')-2s'a-s^S+{s'-5+^/>J)s'_ cosB 



s'(s'-\' X')-2s'cr-s'_5+{s'-5-\' \')s'_ cosB 



- 1 ■ 

-Ui. 



R 



2s' ^( s_;(s'+v^)-2sV-si5-(s'-5+A/y)si cose' 

V^(4 - si C0S6') \s'^{s'-^)-2s'a-s'_5-{s'-S-\f>:)s'_cose ^ 

■ 1 \i _ 2s' 
2-1 i?/ VX! Sdi 



■ I. 



ul 



10.) 



■ 1 ■ 

-U2- 



R 




(s' + V A')-2s'cr+si(5+{s'+<5+V A')si cos6» 



■ lu2 



11. 



12. 



cos6' 



R 



R 



COS 6 



R 



R 



16 s'% 
X' ■ Sdi 



16 s'% 

X' ■ Sd2^ 




^As used throughout the thesis, s (s') is the (reduced) center of mass energy squared, and S12, 
S34 are the invariant masses of final state fermion pairs. The definitions of the derived quantities 
M, ti, i2, Ul, U2, ^1, 22, zi, ^2, si, s^. A, A , A, CT, 5, a* , a| , o" , 02 , 61 , 62 , and c are given in appendices ^ and^ 
where details of the applied phase space parametrizations and the cross-section computation are worked 
out. 

■I'The symmetry operations S and T are defined by: 



S 



/(COS0) 



/(-COS0) 



T g(si2,S34) = ff(s34,Sl2) 
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su ■ (^s — 5 + — sL ■ cos9j (^s — 5 — — sL ■ cos9^ 

As' 

4s[ 

S34 • (s + + VT + si • 0086*) (s + (5 - ^/X + si • cos6') 

As' 

S34 • d2 ^2 



1 



2Si 



H12 



tlt2 



^ ■ + 5 + s'_ cos9^ -S-s'_ cos9^ 



is 
^X 



qt qt 
, ' '-'12+ '-'12- 



ai_ + 6i_ • cos^ 

ai+ + 61+ • cos^ 

02- + &2- • cos^ 

a2+ + 62+ ■ cos6' 



sV - sS^ + -s'__ {s' + 5- 2a) - VX (sS - -s'_ 

-{Vx + s' + s) V 

sV - sS^ + y_ {s' + 5- 2a) + ^/x(^sS- |si) 
+ (v^ -s'-6) -s'_ 

- s^^ + |si {s'-6- 2a) -Vx(^s5+ |si) 

-(yx-s'+5) y_ 

s^a - sS^ + y_ {s'-S- 2a) + ^/X (sS + ^si 

+ (yx + s'-s) y_ 
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.1*1^2- 



R 



■'U1U2 



A1_ 



Au 

AU 
^2+ 



= s 



.tlt2 



Rj 



2S, 



•Inl 



«1U2 



' AU AU 

^1- ^24 
,^1+ ^2- 



2S, 



■I 



UlU2 



8s' 



\f\ + S-s'_ cose^ -S + sL cose^ 



*-'l2+ '-'12- 



8s' 

ai_ — 6i_ • 0086* 

ai+ — 6i+ • cos^ 

02- — &2- • COS^ 

0'2+ — ^2+ • COs6' 



-^1^1- 



R 



^ut 



af 



+ 



R 



1 



= - {s'j^ — s'_ cos9^ — a 



-t2U2- 



R 



af 



af 



+ 



R 



1 

U2 



Rj 



(s| + si cos6l) - (7 = ^(l + 6-cos6l) 



-tiU2- 



R 



12 



In 



2 In 



2A/a 



12 



At 

B2 

Bt 



• hd ■ {a\ + 61) + C12 + 5* ■ (a,i - 6^) /2 
vC^ ■ hd ■ {a\ - h) + C12 + B\ ■ (a,i + bd) /2 

• • (a^ - ^^2) + + B^ ■ (a,i - bd) /2 

• &d • (a^ + &2) + C12 + B^ ■ (a,i + bd) /2 



12 



^12 
-B12 



(Ai2 • COS^^ + B12) 



V 2s' , 
-s'J ■ S12 ■ (s - S12) 

s'J • s • (s — S12) — 2ssi (s — (j) (s — S12) + s^ (s — a)^ 
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Bi 



TDU 

-"2 



-2 
-2 



hi {a\b2 + a^bi) + a^ic^ 



a,i = al + a^ = (s' + 5) - 2s 







= 62- 


" 1 ■ 












-t2Ul. 


R 





2s' 

-27'- 



12 



t\U2 



2VC: 



21 



In I 



■ i?r 



21nf"" + ^^ 



At 

Bt 



a 



21 



B34 



-"2 



as2 



VC^i ■ bd ■ {ai - 62) + C21 + B\ ■ (a,2 - ba) /2 
VCTi ■ bd ■ (a* + 62) + C21 + Bl ■ (a,2 + ba) /2 
■ hd ■ « + 61) + C21 + ■ (a,2 - 6d) /2 
\/^ • bd ■ {al - 61) + C21 + 5" • (a,2 + 6d) /2 

• (/I34 • cos^^ + 534) 



V 2s' , 

-S'J • S34 • (s - S34) 

s'^ ■ s ■ {s — S34) — 2ss'_ {s — a) {s — S34) + {s 
-2 [62 (ai&2 + albi) + 0,20^^ 

-2 bl (0^62 + 461) + as2C^ 

4 + = ^ (S' - - 2S34 



COS 9r 



tl 



R 



16 g'^ bl 
XSdi^ 





■1- 




'1' 


( 




- a{ 

R 


A. 



J R/ 



COS 



^^2 



16 s'2 62 

XSd2^ 



R 
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tit2 



J R 



+ (&2 {a\b2 - alh) - a^c^) 



R 



Ut 



l''2 



03 — 0* = ~2~' ^'^ ~ ^^^^-^ 



^1^2 J 



^ ( 4 - 4(4ad-^2&, 

1*2 V 

- (bi {a\b2 - a^bi) - a\c^ 



R 



Ut 



Rj 



t\t2Ui 



R 



„ut 
"1 



lf'2 



+ 



R 



t2Ui 



RJ 



tit2U2 



R 



af 



tlt2 



+ 



R 



tiU2 



RJ 



1 ■ 






■ 1 ■ 




■ 1 ■ 


tiUiU2. 


R 


af V 


.U1U2. 


+ 

R 


.tlU2. 


1 ■ 






■ 1 ■ 


+ 

R 


■ 1 ■ 


t2UiU2. 


R 




.U1U2. 


.^21*1- 



RJ 



RJ 



1 



tit2UiU2 



R 



1 1 

"1 "2 





■ 1 ■ 




■ 1 ■ 




■ 1 ■ 






-^1^2- 


+ 




+ 


.tiU2. 


+ 






.U1U2. 


R 


R 
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G.4 Bremsstrahlung Integrals — Second Series 



Below, the integrals needed for the last analytical integration in the bremsstrahlung case 
are given. This integration is over the photon scattering angle 9. Using the electron mass 
me, the ultrarelativistic approximation!], and the notation 



[A]g = ]^ j dco^d A , 

-1 



(G.4) 



the following integrals are required for the "second series" of bremsstrahlung integrals: 

1-) 11). = 1 



2 ^ 



3.1 



1 

Z2 



4.1 



s 
4 



5.1 



Z2 



s 
4 



6.1 



S. 



dl 



2s'_ 



7.1 



S, 



d2 



2s'_ 



Ass' 



9.1 



Sd2^ 



Ass' 



10. 



-Zl 

Lci 

Lc5 



D 



zltl 



h -Lti - D 



zltl 



Lc5 ~ Lcl 

s'Js' + v^) + s'a-s6^ 



In 



In I 



Li2 



sL{s' -\^) + s'a- s5 

' s'-a + VX \ 
^s'-a-VYj 

s'_{s + sL - 6 - y/X) 
2s{sL -6) + s'{a + 6) 



Lis 



'Us + s'_-5 + VX') 
2s{sL -6) + s'{a + 6) 



^To obtain the correct result, me rnust be retained in all integrals containing factors 1/zi. In addition, 
keeping me is sometimes suggested for the sake of stability of the numerical invariant mass integrations 
that follow the 6 integration described in this section. 
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[id 

DzU2 
Dz2tl 

1 

— H2 
LZ2 }0 

Lt2 
Lc2 

Dz2t2 

\->A 

— Iu2 

lZ\ J Q 

— lul 

1Z2 J e 

— Iu2 
IZ2 J 



h ■ Lc5 — D 



zlt2 



' -s'_{s' + 6 - ^/X) 
s'{a + S) 



- Lio 



s'{a + S) 



2 ( ■ -^c5 ~ Dz2tl 



-s'_{s' - S - VX) 
s'{a - 6) 



Lis 



-s'_{s' -6 + VX') 
s'{a - 6) 



h ■Lt2 - D 



z2t2 



Lc5 ~ Lc2 



In 



s'_{s' -VX) + s'a + sS 



( sL{s + s'_ + 5-^/X) 
[ 2s{sL + S) + s'{a - S) 



Lis 



s'_{s + s'_ + 5 + VX) 
2s{s'_ + S) + s'{a - S) 
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lZ2 J 



















IZ2 J 













\-'A 
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IZi J 













\-'A 




e 


IZi J 
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1 

Sdi 



tl 



1 

2Z 





\ a + d 



2 s: 



19.) 



S, 



H2 



d2 



1 

2Z 



T 



1 

Sdi 



a 



^. (sis'-S + VX) 





s'{a - 5) 



s'-6-VX' 
a-6 



2 si 



20.) 


iJdl J 


e 




e 


h2 

lJd2 i 


21.) 


\^ 1 ] 


e 




e 


7;— ki 

IJdl i 



22.1 



H12 



h ■ Lci + Dziti + D 



zlt2 



23.) 



— kl2 

Z2 



A comment on the treatment of arguments of logarithms and dilogarithms is in order here. The 
logarithm is undefined for negative real values, the Dilogarithm has a cut for real numbers larger than 
1 as is explained in appendix |^. This means that, in principle, infinitesimal imaginary parts from 
the propagator denominators have to be carried along in the whole calculation. This, however, is an 
unnecessary nuisance, because all integrands are real and regular. If infinitesimal imaginary parts were 
needed, because logarithms or Dilogarithms lay on cuts in final results, they were attributed to the 
invariant masses s\i and S34. As integrands are real and regular, this is a correct treatment, because it 
is then irrelevant for the integrand how an infinitesimal imaginary part is entered. This technique will 
be used in subsequent integrals without further notice. 
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ut 



Hl2 



^c8 



^c6 
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In 



Lc8 ■ -f'ce + -0*1 



s — a 
s' — a , 

(ai_ 6 + 6i_)(a2+ & + 62h 
(ai+ h + (a2- 6 + 62- 



see 15.) & 17.) in app. ^ 



-Li2| 

-Li2| 
+ Li2| 
+ Li2| 
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02+ 
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1 



^c7 



-'-^a2 



In 
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+ Lis 
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-62_(1 
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Li2(^ 
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4s'_s'si2 



'X In \ —L^ 
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■icb 



+ 2^ 



Le3 = In 1 + 



s'_{s-5y 



'd2 



"2 



1 

4sis'S34 



A In — + — — Lc4 — -i^c5 



^c4 



In 1 + 



s'Si2 , 



'dl 



"3 



16 si 



+ 



2 s'^Si2 



^c3 



+ 



(7^ - (5^ + 2(7(^ + 2s'{a -S)- 2s''' 



^c5 



„/2 



"27^^2 



12 



'(i2 



16 sL 



+ 
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SS'^S34 2 s'^S34 



S'^S34 



In h 

s 



2 S'^S34 
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40.) 



cos^ 



1 ( s-S 



^cl 



41.) 



cos^ 
dt d2 



1 (s + 5 



^c2 



42.) 



ct ct 

.'-'12+ '-'12- 



4si\/A 



= In 



+ f7 + yA^ 

+ (7- V^, 



43.) 



cos^ 



4s: 



/ 2 



, S34 VA 
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- Li2 
+ Li2 

- Li2 
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- Li2 

+ Li2 



c_| a e+ 
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c_|__a+_e_ 
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+ Li2H 
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s' -a±VX 

8 S Si2 S34 



- Li2 
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d 


c_ 





103 



H12 



Qt Qt 2 
.'-'12+ '-'12- 
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^tl2 


J. " 2A 


qt qt 

.'-'12+ '-'12-. 



+ 



8sLX 



X'VX ^ - ^ r s'_ + S ^ 
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2SS12S34 ^ 2 S^Si2 

A - s(s!_ + a) 



2 sis34 



4SS12S34 

4sSi2S34 



(Lci -h Lc2) 



(-f'c3 + Lca) 



C0S9 ■ ltl2 


S 


^tl2 


qt 2 qt 2 
.O12+ '-'12- . 


e ~ 


qt "2 qt 2 
.O12+ '-'12- . 



+ 



-^cl + TT"? -^c2 



2s'Si2S34 S34 2sisi2 2sis34 



S — (T 

4 SS12S34 
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4 SS12S34 
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^tl2 
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^tl2 
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_Oi2+ '-'12- . 




ct qt 2 

_Oi2+ '-'12- . 
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32siv^' 



■X 
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A — '1 — 9 — 9 ~A — 9 — 9~ -f'ci2 + — ; — 9~ L^i 

4ssis|2S34 4sf2S34 4slsf2 



^/X VT (2SS12S34 + A(s' - cr 



SS12S34 - + (7) + A(s| - a] 

8s^Sl2'534 

\/A vT" (s(si + a) -2 S12S34 - A 



+ 



{Lcl + Lc2) 



8s^Si2S34 



(Lc3 + Lci) 



^cl2 — 



1 (si si, 5^/X\ 

— ■ \ — ^c2 — r 

si y si SI s J 
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50.1 



cosO ■ I 



tl2 



at Q 
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4A 



COS^ ■ Itu 
at 2 2 

^12+ '-'12-- 
^tl2 



Qt 3 Qj 3 
L^12+ '-'12- J 



+ 



4siA 



+ 



32 sL'X 
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H12 
Qt 2 Qt 2 
L^12-|- '-'12- 



^2 „2 -^cl2 



+ 
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2 si si2 ''^ 2 si S34 
(^(s'-cr) A 2SS12S34 

8 S^sf2S34 

f21n£H 

V S34 

AA' + 2 S12S34 (s (s'-cr) + A 

8 e2,2 „2 
OS S]_2'534 



c2 
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Lc3 + Lei 



cl — ^c2 



51.) 



''tlU2 



2v^ 
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1 • C12 
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/3i2 + ixo \ 
(3i2 + ia;o 



+ Li2(^ --,— j + Li2 



Li2 
Li2 



' /?i2 + ixp ' 
. ^2 . 




- Li2| 
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C12 



/9: 
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Tl 



T2 
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(ss' - SS34 - S'si2' 



S_ S12 ■ Xq 



V\' S S12 ■ Xo 

S12 (s - S) 



^^Integrals 51.) to 60.) were not directly calculated from the result of the il/j integration. Instead, a 
Feynman parametrization was used to "linearize" l/tiU2 and 1/^21*1: ^/tiU2 = Jq da/ [tia + U2{1 — a)]^- 
Then, Qn, 9, and finally the Feynman parameter a were integrated. 
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fc+i 
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> 
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^Hcrc and in the remainder of this section the expressions obtained by the symmetry operation T will 
not be rewritten. The reader can easily do that by making the reciprocal replacement S12 S34. 
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/° 

''1 




''2 
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= T2 


^4 


= — T2 


^5 


— T* 


''6 


— — T* 

— ^2 


h 


T-* 

- ^1 




T-* 

- -^1 


^k± 








At) 






_ ia;o + 



Xl/2 = 



si Sl2(s - 5) ± \fX {ss' - SS34 - s'si2) 

72 V 
si S12 — sX 



Ami 



54.) 



21 



tlU2 



• 2 \/C\2 



4 (SS' - SSyi - S'S34) 



<^ 2 • (2Lc5 — 



55.) 



''t\U-2. 



Z2-2y/C: 



12 



^i-2V^ 



12 



56. 



''t2'f 1 



^2-2V^ 
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57.) 



12 



— ■ D" 

o 



For the calculation of the integrals Xi the following quantities are needed: 



A3 
S3 

A. = A'- 



si2(s;-2a)2 - ^g^sA' 
S12 (s+ - 2(t) (s - 5) 

|ss' - SS34 - s'si2 - 2si2 (s - a)| 

- Ami {si2 (s - + a' (s - S12)} 
-53±x/A; 



In principle, depending on the values of 2:3 and 0:4, four cases have to be distin- 
guished: 

(i) Xs/^eM; {x^l > {x^l > [3 
(ii) X3/4^e]R; |x4| > P > \x3\ 
{in) 2:3/4 eiR; P > \x4\ > \xs\ 
(iv) A3 < 0; xs = xl e € 

However, since the integrand is regular, it is sufficient to present the solution for 
case (i). The solutions for the other three cases are then obtained by analytical 
continuation. In case (iv). which is relevant for only a very small fraction of 
the phase space, the problem of a numerically correct treatment arises. This 
problem was solved by computing the integral's value for a very nearby phase 
space point so that the expression of case (i) could be used. Having in mind that 
the integral is regular, it is clear that thus only a negligible error is introduced. 
The expressions Xq, Xi, and X2 for case (i) are presented below. 



Xn 



^c8 



In I 



- Inl 



' P — X4 

+ X4^ 
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"\/si2 • (s[ - cr) 
2 VsA' As~ 



E(-i 

k=l 



(-1)^+^ • i^i - «34) 
i=3 yjxf — Xq 

Ai) 



In 



+ Li. 




+ Lin 



Lio 




'4 - 2(7 

-(0 



X7, 



ixo + ^x^ — x^ 



t{x) - 



X 



2x/A^ 

2 2 



X 



i=i 3=1 



In 



~ ' R 

X j — /3 



37 ^ 37 j 1£ 



+ Li2 



X X j 1£ 



+ Li2 



Xj^ I ^'"^ 



- Lio 



37^ I 
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G.5 Loop Integrals 



For the computation of the virtual initial state QED corrections to the studied pro- 
cess apart from the integration over final state fermion decay angles, two sets of 
integrals are needed. The second set is given in section |G.(j| . The first set is presented in 
this section and consists of the integration over the loop momentum p. 

It is sufficient to present the integrals for the virtual t-channel graphs, because the in- 
terferences of the virtual u-channel graphs with the Born level t- and u-channel graphs 
equal the interferences of the virtual t-channel graphs with the Born level u-channel and 
t-channel graphs. This is easily verified by symmetry arguments. First, some definitions 
are introduced: 



1 

1 


1 

p"^ — ie 

1 




1 


1 


{p — kiY + ml 
1 


— ie 


p2 — 2pki — ie 
1 


1 


{p + + ml 
1 


— ie 


p2 + 2pk2 — ie 




(p - QtY + ml 


— ie 





(G.5) 



The momenta ki, /c2, and qt are carried by the incoming electron, the incoming positron 
and the intermediate initial state electron. For brevity, q = qt with q^ = t (compare equa- 
tion (|E.5|) ) will be used subsequently. Dimensional regularization is used to consistently 
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treat divergent loop integrals. The following notations are employed: 



and 



(27r)'* 



A 



IJ.=me 



n — 4: 2 II ^ ' 



with |n - 4| < 1 , 
n = 4 — £ , 
n = 4 + £ , 



IJ-=me 



tJ-=me 



= ln| 
^2 = ln| 



— 1£^ 



In — - — 'm = la — in 



Sl2 



1 Sl2 . 

9 • ; = m— - ITT 

mj — IS J TTlg 



S34 



— i£ 



1 S34 

m — - — ITT 

ml 



(G.6) 
(G.7) 
(G.8) 



(G.9) 



Euler's constant has the value je = 0.577216... . The difference between the ultraviolet 
divergence P and the infrared divergence P^^ is that P is regularized by n = 4 — £ and 
P^^ by n = 4 + £ with < £ < 1 . In the URA one obtains: 



1-) 



3.1 



noni 



.noni 



i-fcf 

167r2 



P + 



2.) 
4.) 



M0II2. 



167r2 



P + 



5. ) 

6. ) 

7. ) 

8. ) 



non„ 



• 1 



= — ^ ( -2P-lt + 2] 



167r2 



^(-2P-Zo + 2) 
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11.1 



16 tt' 



12. ) 

13. ) 



-noHina 



_ ^ piR + i f ^ _ 

167r2 s ' s I 6 2 



14. 



-noHina 



16 TT 
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167r2 



^ \ 2 A A) 
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16.) 
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1 


167r2 


t + S34 
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noHin, 
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P 3 1 
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t 



2 t + S12 
3^2 
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t + S12 



2(i+Si2) \t + Si2 



{k - h) 

1 

0* - /i) - 1 
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2(t + Si2) 



(^t - ^l) 



non2n. 
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t + S34 
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(t ■lt + S34 ■ I 



1 

2 t + S34 
3^2 



{t + 534)- 
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16 TT^ 



X1-X2 



2s s 

1 — X2 s + S + y/X 
1-Xi s + S - VX 



167r2-A 



X 



(s - a) Ii2q - {s + d)ln^ - (s-5)ln^j 
+ /cf ^ - S34 (^s + Ii2q + (s - (j) In^ + 2 S34 In^ 
+ k2 ^Si2 (s - ^) Ii2q - 2 S12 - (s - (j) lii-^ j 
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26 



1 
A 
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SS34 



s + S34 6 SS12S34 
2 A 

■^34 3 SS34 

~2 A~ 

3Sl2 

A 



■In 



S12 
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■I 



12q 
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+ 



■In 
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As mentioned in appendix Pj^.2| , a FORM computer algebra program was used to handle 
the many terms in the virtual corrections' calculation. Inspection of the FORM code 
at an intermediate level showed that, for the calculation of virtual QED corrections to 
process (|1 . 1| ) (see appendix [E.2|) , the four-point function is only needed in one specific 
Lorentz contraction. Precisely this Lorentz contraction is given below. 
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nonin2n„ 
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k {h + h 



21 
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G.6 Virtual Corrections' Phase Space Integrals 



After the integration over the loop momentum, there is one non-trivial integration left 
when the virtual initial state corrections to process ( [I.ID are to be evaluated. This is 
the integration over the boson scattering angle i) in the center of mass system. It should 
be recalled that it is sufficient to present the integrals for the t-channel case and the 
t-channel / u-channel interference for reasons of symmetry. Making use of the URA and 
the notations 
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the following list of integrals is obtained. 
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